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Abstract 

Analytical models for IEEE 802.11 -based WLANs are invariably based on approximations, such as the well-known 
decoupling approximation proposed by Bianchi for modeling single cell WLANs consisting of saturated nodes. In 
this paper, we provide a new approach to model the situation when the nodes are not saturated. We study a State 
Dependent Attempt Rate (SDAR) approximation to model M queues (one queue per node) served by the CSMA/CA 
protocol as standardized in the IEEE 802. 1 1 DCF MAC protocol. The approximation is that, when n of the M queues 
are non-empty, the transmission attempt probability of the n non-empty nodes is given by the long-term transmission 
attempt probability of n "saturated" nodes as provided by Bianchi's model. The SDAR approximation reduces a single 
cell WLAN with non-saturated nodes to a "coupled queue system". When packets arrive to the M queues according 
to independent Poisson processes, we provide a Markov model for the coupled queue system with SDAR service. The 
main contribution of this paper is to provide an analysis of the coupled queue process by studying a lower dimensional 
process, and by introducing a certain conditional independence approximation. We show that the SDAR model of 
contention provides an accurate model for the DCF MAC protocol in single cells, and report the simulation speed-ups 
thus obtained by our model-based simulation. 

Keywords: IEEE 802. 1 1 DCF, CSMA/CA, single cell, coupled queue system, state space reduction, iterative 
method, model-based simulation 



1. Introduction 

The IEEE 802.11 standard J2l has been widely adopted as the de facto standard for accessing shared wireless 
media in Wireless Local Area Networks (WLANs). The Distributed Coordination Function (DCF) Medium Access 
Control (MAC) protocol, which is a particular version of the Carrier Sense Multiple Access with Collision Avoidance 
(CSMA/CA) protocols, provides the fundamental access method in 802.11 WLANs D2J]. The Hybrid Coordination 
Function (HCF) which provides service differentiation and Quality of Service (QoS), and the optional Point Coordi- 
nation Function (PCF) are also built on top of the DCF |2J|. Hence, understanding the behavior of the DCF is the key 
to understand the performance of 802. 1 1 based WLANs. 

This paper is concerned with analytical modeling of DCF-based WLANs in the single cell scenario. As in i20ll . 
we define a single cell to be a set of closely located 802. 1 1 nodes such that: (i) every node can sense the transmissions 
by every other node, and (ii) every node can decode the transmissions by every other node in absence of interference. 
Figure Q] depicts an example of a single cell WLAN which operates in the ad hoc mode. In this example, users or 
client stations (STAs) can directly communicate among themselves. Figure |2]depicts another example of a single cell 
WLAN which operates in the infrastructure mode. In this latter example, STAs are primarily interested in accessing 
the Internet through an Access Point (AP). The analytical model in this paper applies to single cells that operate either 
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Figure 1: A single cell in the ad hoc mode: In this 
case, users or client stations (STAs) can directly com- 
municate among themselves. 



Figure 2: A single cell in the infrastructure mode: In this case, STAs commu- 
nicate among themselves and access the Internet through an Access Point (AP). 



in the ad hoc mode or in the infrastructure mode. To develop our analytical model, we do not need to distinguish 
between the AP and the STAs, and we call an AP or a STA, a "node". We assume that the readers are familiar with 
the DCF MAC protocol, the details of which can be found in the standard document |2|]. 

Analytical modeling of 802.1 1 WLANs has been a topic of great interest ever since the publication of the original 
version of the standard OJJ] . In his seminal work, Bianchi [5] proposed an accurate analytical model for DCF-based 
single cell WLANs consisting of saturated nodes. A node is said to be saturated if it always has a packet to transmit, 
i.e., if its transmission queue never becomes empty. The saturation model is often useful in providing a sufficient 
condition for stability |3, 22] Q In fact, the saturation assumption adopted by Bianchi [5] was partly motivated by 
the requirement to determine the maximum traffic intensity that the system can sustain under stable conditions. The 
saturation assumption, of course, simplifies the modeling problem since the node queues never become empty and 
queueing dynamics can be ignored. However, the saturation assumption is not valid, in general. In many real applica- 
tions, e.g., web browsing and packetized voice telephony, a single cell WLAN might operate well below saturation. 

When a node is not saturated, its (transmission) queue can be empty for a positive fraction of time during which 
it does not contend for accessing the medium. Thus, queueing dynamics must be studied in order to determine the 
fraction of time for which a node contends. Due to the contention for medium access, however, the queue length 
processes of the nodes are coupled. Hence, modeling a single cell WLAN with non-saturated nodes involves ana- 
lyzing a "coupled queue system " where the queues are served according to the CSMA/CA protocol as standardized 
in IEEE 802.11 DCF yfl. In this paper, we apply a State Dependent Attempt Rate (SDAR) approximation (see 
Approximation lXTl in Section [3T2l > to study such a coupled queue system. 

1.1. Literature Survey 

Bianchi ^ developed the seminal analytical model for DCF-based single cell WLANs consisting of saturated 
nodes. The key approximation in Bianchi's model is the following: 

Approximation 1.1 (Bianchi's Approximation). Every transmission attempt by a node collides with a constant prob- 
ability 7 regardless of the history of collisions. ■ 

Approximation [I"T1 is popularly known as the decoupling approximation since it decouples the backoff processes 
of the nodes through a constant parameter 7. In reality, the backoff processes of the nodes are coupled precisely due 



'it is shown in Q] that, if the arrival point process is a stationary and ergodic marked point process of finite intensity, then under certain 
constraints (which are satisfied by several classical open queueing networks), the system is stable if the intensity of the arrival process is strictly 
smaller than the "intensity of the departure process obtained by saturating the queues." By "stability" here we mean that the (joint) queue length 
process has a proper stationary distribution. 
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to collisions. Applying the decoupling approximation, Approximation ITTTl Bianchi considered a tagged node, and 
developed a two-dimensional Markov chain to model the evolution of the backoff process of the tagged node. Since 
the backoff parameters (i.e., CW m i n , CW max , and the retransmit limit) of the nodes are identical in the DCF, any 
node can be taken as the tagged node. The two-dimensional Markov chain keeps track of the backoff stage (i.e., the 
number of times the same packet has already collided) and the backoff count (i.e., the remaining number of backoff 
slots for the next transmission attempt to begin) of the tagged node. From the stationary distribution of the two- 
dimensional Markov chain, Bianchi obtained "the probability (3 that a node attempts a transmission in a randomly 
chosen slot" as a function of the collision probability 7. Let 

(3 = G( 7 ) 

denote this function. With n saturated nodes in the single cell, Bianchi obtained the collision probability 7 as a 
function of the attempt probability j3 by 

7 = r(/3) : =1-(1 -/J)"-*, 

since every node attempts a transmission in any randomly chosen slot with probability j3. We emphasize that, 7 is not 
the probability that a collision occurs in a randomly chosen slot; 7 is the probability that an attempted transmission 
collides. The above two equations yield a fixed point equation 7 = T(G(^)). The transmission attempt probability j3 
and the system throughput could be obtained using the solution of the fixed point equation. 

Cali et al. 10, Ht] proposed and analyzed a p-persistent version of IEEE 802.11. Considering saturated nodes, 
they obtained a theoretical upper limit, and also proposed an adaptive backoff mechanism to achieve the theoretical 
throughput limit. Kumar et al. [20] generalized the Bianchi model to arbitrary distributions of backoff times, backoff 
multipliers and retransmit limits. 

Motivated by the need to understand the performance of WLANs under realistic traffic conditions, modeling of 
the non-saturated case has attracted much attention in the recent past. Several models for the traffic arrival processes 
have been considered in the literature. Under the "Poisson traffic" model, packets of fixed size are assumed to arrive 
into the node queues according to independent Poisson pro cesses of given rates oil I9I EM [lol fl, [l6l 14, 13\. Models 



for packetized voice telephony have been considered in jl5tl23ll . In OA the authors consider ON-OFF traffic sources 
with exponentially distributed OFF periods, geometrically distributed number of packets during ON periods, and 
exponentially distributed packet payload sizes. The Poisson, voice, and ON-OFF traffic models belong to the so- 
called open-loop type where the traffic sources do not adjust their sending rates depending on the level of congestion 
in the network. Arrivals according to the closed-loop control of TCP have been considered for long-lived flows 
in and for short-lived flows in 1271 12411 . 



As pointed out earlier, modeling the non-saturated case requires analyzing a coupled queue system. To that end, 
several simplifying assumptions have been made in the literature, either explicitly or implicitly. A detailed account 
of the veracity of common modeling hypotheses, in the context of 802.1 1 WLANs, can be found in lfl7ll . We follow 
the precise terminology of 11711 to comment on the relevant assumptions. For any tagged node, define Ck '■= 1 if the 
k th transmission attempt by the node results in a collision, and define Ck '■= otherwise. Also, for any tagged node, 
define Qk := 1 if there is at least one packet awaiting transmission after the k th successful transmission from the 
node, and define Qk := otherwise. The simplifying assumptions in 113 lfl — 1 1 Oil are then equivalent to the following 
assumptions regarding the sequences {Ck} and {Qk}' 

(Al) The sequence {Ck} consists of independent random variables. 

(A2) The sequence {Ck} consists of identically distributed random variables. 

(A3) The sequence {Qk} consists of independent random variables. 

(A4) The sequence {Qk} consists of identically distributed random variables. 

In particular, the authors in ll3lll - lfloll assume that: (i) in a randomly chosen slot, each node attempts a transmission 
with constant probability f3, (ii) every (transmission) attempt collides with constant probability 7, and (iii) each queue 
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is not empty after a departure with constant probability They analyze the evolution of the backoff process and the 
queue length process of each node in isolation and obtain fixed point equations relating f3, 7, and q. The probabilities 
(3, 7, and q are obtained by solving the fixed point equations. The collision probabilities, the throughputs and the 
mean packet delays are obtained using the solution of the fixed point equations. 

Note that, Assumptions (A1)-(A4), essentially "decouple" the transmission attempt processes, the backoff pro- 
cesses and the queue length processes of the nodes through the (unknown) constant probabilities (3, 7, and q, which are 
obtained by solving certain fixed point equations. Assumptions (Al) and (A2) regarding the collision sequence {C&} 
constitute the decoupling approximation introduced by Bianchi which can be called a collision-decoupling ap- 
proximation. Assumptions (A3) and (A4) regarding the queue-status sequence {Qk} are specific to the non-saturated 
case and, together, they can be called a queue-decoupling approximation. 

In |0], the authors assume only (Al) and (A2) and apply a matrix-geometric analytic method. It was shown 
in lfl6ll that Assumption (A4) leads to inaccurate predictions for throughputs when the arrival rates into the queues 



differ significantly from each other. In [14] it is argued, by providing results from NS-2 simulations 02611 . that the 



collision-decoupling approximation of Bianchi works well in the saturated case, but leads to inaccurate results in the 



non-saturated case. A detailed discussion of the validity of Assumptions (A1)-(A4) can be found in [17] where the 
authors provide evidence from NS-2 simulations to conclude that: 

1. Assumption (Al) is valid regardless of whether the nodes are saturated or not. 

2. Assumption (A2) is valid when the nodes are saturated but is not valid, in general, when the nodes are not 
saturated. 

3. Assumption (A3) and (A4) are not valid. 

In 014ll . the authors propose to model the transmission attempt probability of the nodes at any instant t as a function 



of the number of non-empty nodes in the system at t. Clearly, in 01411 . the attempt probabilities of the nodes are state- 
dependent. The approach of state-dependent attempt probabilities has also been adopted in ll23l l6Tl. In 1 1 3 , 27 , 24 1 a 
state-dependent service rate approach is adopted. 

1.2. Preview of Contributions 

We develop a new approach to model single cells with non-saturated nodes under Poisson packet arrivals. Guided 
by the reported inaccuracy of the state-independent approach, we adopt a state-dependent attempt probability approach 
as in ll 14i I23L I6TI . The state-dependent attempt probabilities in lfl411 are obtained by an iterative method which requires 
computations involving a three-dimensional Markov chain. We, however, apply an approximation proposed in B23I1 to 
obtain the state-dependent attempt probabilities (see Approximation [3~T1 in Section EOl i. As explained in Section [5721 
our model is computationally less expensive than that in I1 1411 . In particular, our model requires computations involving 
a two-dimensional Markov chain. We emphasize that, even though we apply the approximation proposed in B23I1 . the 
problem we address in this paper is completely different from that in 02311 . The problem setting in 02311 is such that 
analysis of queueing dynamics is not required whereas we analyze a coupled queue system. 

Our contributions in this paper are the following: 



We develop a Markov model with Poisson packet arrivals. Our Markov model reduces a single cell WLAN 
with non-saturated nodes to a coupled queue system with SDAR service discipline (Section |5J. We provide a 
sufficient condition under which the joint queue length Markov chain is positive recurrent (Theorem l3.U . 



For the case when the arrival rates into the queues are equal, we propose a technique to reduce the state space of 
the coupled queue system (Section|4]i. For the case when the buffer sizes of the queues are finite and equal, we 
propose an iterative method to obtain the stationary distribution of the reduced state process (Section l57fT i. Our 



iterative method is computationally less expensive than that in [ 14] (Section [5T2J, and yet, it provides accurate 
predictions for important performance measures (Section|7j. 



2 In general, if the nodes are not identically parametrized, e.g., in 802. lie, then Node-i is associated with the probabilities ft, 7$, and qi. 
However, for a given node i, the probabilities ft, 7;, and qi are assumed to be constants independent of the current state of the system. 
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• We applied the SDAR approximation to modify the MAC layer of NS-2, keeping all other layers unchanged. 
Originally, our objective in doing so was to validate the SDAR approximation itself by comparing the results 
obtained from the unmodified and modified NS-2 simulations. However, by doing so, we could also demonstrate 
the possibility of improving the speed of simulations by model-based simulation at the MAC layer. We show 
that the SDAR model of contention provides an accurate model for the CSMA/CA protocol in single cells and, 
at the same time, achieves speed-ups (w.r.t. MAC layer operations) up to 1.55 to 5.4 depending on the arrival 
rates and the number of nodes in the single cell WLAN. 

1.3. Outline of the Paper 

We summarize our network model and assumptions in Section [2] In Section [3] we introduce the SDAR approxi- 
mation and develop a Markov model with Poisson packet arrivals and infinite buffers. In Section |U assuming equal 
arrival rates, we reduce the state space of the coupled queue system and obtain the transition probability matrix of 
the reduced state process. To demonstrate the predictive capability of our model, we restrict to the case of finite and 
equal buffers in Section [5] In Section [57X1 we propose an iterative method to obtain the stationary distribution of the 
reduced state process (for equal arrival rates, finite and equal buffers), using which, we obtain predictions for impor- 
tant performance measures in Section l5~3l In Section|6]we report how the SDAR heuristic technique could be applied 
to improve the speed of simulations. In Section [7] we validate our coupled queue model and our iterative method by 
comparing with NS-2 simulations where we also discuss some simulation results for the case of unequal arrival rates. 
Section[8]concludes the paper. Proofs and derivations have been provided in the appendices. 

2. Network Model and Assumptions 

We consider a IEEE 802.11 DCF-based single cell WLAN consisting of M nodes. Our network model and 
assumptions are given in the following: 

• The nodes are homogeneous. This means that the nodes use identical backoff parameters, i.e., they use identical 
CWmtn, CW max , and "retransmit limit" (which is indeed true for the DCF MAC protocol). 

• The arrival processes bring packets of fixed size L into the node queues according to independent Poisson 
processes. The arrival rate in packets/sec into Node-z's queue is denoted by 

• The nodes use equal Physical layer (PHY) rates to transmit their packets. 

• The wireless channel is error-free which implies that single transmissions are always successful. 

• There is no packet capture. Simultaneous transmissions (i.e., collisions) always result in the failure of all the 
involved transmissions. Thus, there can be at most one successful transmission at any point of time. 

3. A Coupled Queue Model 

Since there can be at most one successful transmission at any point of time, the system can be viewed as a single 
server serving multiple queues. The M queues corresponding to the M nodes are coupled essentially due to MAC 
contention and our immediate objective is to model the evolution of the joint queue length process. We proceed by 
embedding at the so-called "channel slot boundaries" described in the following. 

3.1. Channel Slots 

Let a denote the duration of a backoff slot in seconds. The duration a is a PHY parameter Q]. We call a time unit 
equal to the duration a of a backoff slot, a system slot. We call a period of time during which all the nodes in the system 
have empty queues, a system-empty period. For analytical convenience, we make the following approximations: 

(ai) Nodes always sample non-zero backoffs. Consequently, when an activity period ends (i.e., after a successful 

transmission or a collision), subsequent transmission attempts can occur only after at least one backoff slot. 
(a-i) System-empty periods are integer multiples of system slots. 
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Figure 3: Description of channel slots: The channel slot boundaries indicated by arrows are the only possible attempt instants. 

As per the rules of the DCF, in reality, nodes do sample backoffs with some positive probability. Approximation 
(0,2) is also not true, in general, since arrivals occur in real time. However, it will be clear from the accuracy of our 
analytical model that the errors due to Approximations (ai) and (02) are negligible. Owing to Approximation (ai), 
transmission attempts in the system can possibly occur only immediately after the end of a backoff slot. These possible 
attempt instants have been indicated by arrows in Figure [3] which depicts the backoffs and the activities in a single 
cell. The channel slots that occur on the common medium have also been shown. We call the time interval between 
any two such possible attempt instants a channel slot, and observe that the channel activity evolves over cycles of 
channel slots. Note that, channel slot boundaries are the only possible attempt instants. 

When the system is empty, an idle channel slot of duration a occurs (see Approximation (02)). A succession of 
idle channel slots occur until arrivals make some of the nodes non-empty. Non-empty nodes sample non-zero backoffs 
(see Approximation (ai)) and attempt transmissions when their backoff counters become 0. Depending on whether 
there are no attempts, only one attempt, or more than one attempt made in the system, an idle, a success, or a collision 
channel slot occurs. The duration of an idle channel slot when the system is non-empty is equal to the duration a of a 
backoff slot. By Approximation (02), when an activity period ends, subsequent transmission attempts can only occur 
after a backoff slot. Hence, we combine the time duration of a seconds, which immediately follows an activity period, 
with the activity period itself to form success or collision channel slots. The attempt process resumes at the end of 
channel slots, thereby creating more channel slots and the process repeats. 

Let T s (resp. T c ) denote the duration of a successful transmission (resp. a collision). The success time T s and 
the collision time T c depend on the PHY rate of transmission. Since the nodes use equal PHY rates and the packets 
are of fixed size L, the durations of T s and T c are fixed and equal for all the nodes. The duration T s corresponds 
to "DATA-SIFS-ACK-DIFS" in the Basic Access mode, or "RTS-SIFS-CTS-SIFS-DATA-SIFS-ACK-DIFS" in the 
RTS/CTS mode. The duration T c corresponds to "DATA-DIFS" in the Basic Access mode, or "RTS-DIFS" in the 
RTS/CTS mode. Note that T s and T c need not be integer multiples of system slots. The duration of a success channel 
slot is equal to T s + a and that of a collision channel slot is equal to T c + a. 

3.2. Coupled Queue Formulation 

We model the evolution of the system over discrete time instants embedded at the channel slot boundaries. Fig- 
ure 2] depicts the evolution of the queue length process of a typical node in the system. Let T(t),t = 0, 1, 2,3,.. ., 
with T(0) = 0, denote the channel slot boundaries. The t th channel slot (t > 1) is precisely the time interval 
[T(t — 1), T(i)). The duration of the t th channel slot is denoted by L(t). We analyze the system in discrete time, 
where the discrete time index t corresponds to the actual (i.e., continuous) time instant T(t). 
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Figure 4: Evolution of Node-i's queue. 



Let Qi(t),t > 0,i = 1,2,..., M, denote the number (of packets) in the i th node's queue at time t. Let Ai(t) 
(resp. Di(t)), t > 1, i = 1, 2, . . . , M, denote the number of arrivals into (resp. departures from) the i th node's queue 
in the t th channel slot. Notice the embedding of Qi (t), (t) and D L (t) in Figure[4] Departures from the queues occur 
at the end of channel slots, since a packet is removed from a transmission queue only when an ACK is successfully 
received or a timeout occurs. Arrivals that occur during a channel slot are taken into account only at the next channel 
slot boundary, since, if a node is empty in the beginning of a channel slot, it can attempt only at the next channel slot 
boundary provided that packets arrive into its queue during the channel slot. The queue lengths are updated after the 
arrivals and the departures in the previous channel slot have been taken into account. 

In this section, we assume that each node has infinite buffer spaced Thus, Vi > 0, i = 1,2,..., M, we have 
Qi(t) G N where N := {0,1,2,...}. Also, Vt > l,i = 1,2,..., M, we have A, (t) £ N, and AO) 6 {0,1}. 
The last constraint follows from the fact that there can be at most one successful transmission, and thus, at most one 
departure, in a channel slot. Clearly, the "number in the queue" processes {Qi(t), t > 0}, 1 < i < M, evolve as: 

Qi(t + 1) = Qi{t) - D,(t + 1) + Ai{t + 1). (1) 

Evidently, it must hold that (Qi(t) = 0) =$> (Di(t + 1) = 0), since there cannot be a departure from an empty queue. 

Due to the "Poisson arrivals" assumption, the distribution of the number of arrivals in a channel slot depends only 
on the duration of the channel slot. The duration of a channel slot is known if the channel slot type (i.e., whether it is 
an idle, a success or a collision channel slot) is known (see Section[3j}- Let Lidie, L succ and L co u denote the duration 
in seconds of an idle, a success and a collision channel slot, respectively. Then, we have Lidie = c, L succ = T s + a 
and L co u = T c + a. For a given PHY, the duration a is known. Also, given the packet payload size L and the PHY 
layer (transmission) rate, the durations T s and T c can be computed f5i l20ll . With slight abuse of notation, we indicate 
the occurrence of an idle, a success and a collision channel slot by L(t) — Lidie, L{t) = L succ and L{t) = L co u, 
respectively, and define, Mt > 0, 1 < i < M, the following probabilities: 

di(j) := P(Ai(t + 1) = j | L(t + 1) = L ld ie) 

Vj < 0. 

Si {j) := P(Ai{t + 1) = j | L(t + 1) = L succ ) 

,->,(r. + ,) (*m + .))' ¥jS0> 

V j < 0. 
d(j) := P(A l (t + 1) =j \ L(t + 1) = L coU ) 



e~^^- Vi>0, (2) 



e -A. t yi c ^a) v - v - >_!_ V?>0, 

Vj < 0. 



3 The infinite buffer assumption will be dropped in Section|3]and we will analyze the finite buffer case as well. 
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As described in Section [3TI nodes can attempt only at the channel slot boundaries. Only those nodes that are 
non-empty at t can attempt a transmission at t. Let N(t) denote the number of non-empty nodes in the system at t. 
Then, by definition, we have 



A I 



N (t) =J2 1 {Qdt)>o}, (5) 

i=i 

where lr.i denotes the indicator function. We now introduce an important approximation regarding the attempt 



processes of the nodes which was first proposed in 12311 . 



Approximation 3.1 (State Dependent Attempt Rate). At any channel slot boundary t, t > 0, every non-empty 
node attempts a transmission with probability /3jv(t) where f3 n is the attempt probability of the nodes in a single cell 
consisting of n homogeneous and saturated nodes. ■ 



We call Approximation 13 . 1 1 the State Dependent Attempt Rate (SDAR) approximation. The /3„'s in the SDAR 
approximation can be obtained by a saturation analysis as in J^lllSl- As a consequence of the SDAR approximation, 
given N(t) = n, the number of transmission attempts made in the system at the channel slot boundary t is binomially 
distributed with parameters n and /3„. Hence, the probability that the (t + l) th channel slot is an idle, a success or a 
collision channel slot can be obtained as follows: 



Pidle,n 
Psucc,n 
Pcoll,n 



P(L(t + 1) = L ldle | N(t) = n) = (1 - fi n ) n 

P(L{t + 1) = L succ | N(t) =n) = n/3 n (l - /3„) n_1 \ (6) 

P{L(t + 1) = L coU I N(t) = Tl) = 1 - Vidlem ~ Psucc,r, 



Furthermore, in case of a success channel slot, the packet departure can occur from any of the non-empty queues with 
equal probability. Thus, for all t > 0, 1 < i < M, the number of departures Di (t + 1) from the i th queue in the 
(t + l) th channel slot satisfy 

P (Di(t + 1) = 1 | N{t) - n, L(t + 1) = L succ , Qi(t) > 0) = - . (7) 
1 n 

Note that, if L(t + 1) ^ L succ or Qi(t) = 0, then D t (t + 1) = 0. 

Owing to the channel slot structure imposed by Approximations (ai) and (02), the "Poisson arrivals" assumption, 
and the SDAR approximation, the joint queue length process {Q(t),t > 0}, where 

Q(t) := {Qi(t),Q 2 (t),...,Q M (t)), 
is an A/-dimensional Discrete Time Markov Chain (DTMC) embedded at the channel slot boundaries. 
Theorem 3.1. The DTMC {Q(t), t > 0} is positive recurrent if, 1 < i < M, we have Xi > 0, and 

< min O sa t,n, 

l<n<M 

where Q sa t,n is the aggregate throughput in packets/sec in a single cell consisting of n homogeneous and saturated 
nodes. 

Proof: See Appendix lAl ■ 

Remarks 3.1. The variation of Q S at.n with n, in general, depends on the backoff parameters. However, for the default 
backoff parameters as prescribed in the 802.11 DCF, we observe that, for M large enough (M > 5 suffices), we 
have, mini<„<M ©sat.n = ®sat,M- Then, the DTMC {Q(t), t > 0} is positive recurrent if the aggregate arrival rate 
J2i=i \ i s strictly less than the aggregate throughput <d S at,M of M saturated nodes. Thus, Theorem B . 1 1 exhibits the 
connection between the saturation throughput and the maximum stable throughput of the system. ■ 
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By Theorem im there exist > 0, 1 < i < M, such that the DTMC {Q(t),t > 0} is positive recurrent. We 
denote the stationary distribution of the DTMC {Q(t),t > 0} by v(k), k = (ki, k 2 , ■ ■ ■ , k M ) £ N M . In principle, 
important performance measures such as collision probability, throughput and mean packet delay can be obtained 
once the stationary distribution of the DTMC {Q(t),t > 0} is known. However, the DTMC {Q(t),t > 0} has a 
state space N M , which is difficult to handle for M > 1. The transition structure of the DTMC {Q(t), t > 0} is such 
that, for M > 1, it is difficult to obtain closed-form expressions for the 2/(fc)'sQ Furthermore, it is not possible to 
numerically solve the infinite number of balance equations of the DTMC {Q(t), t > 0}. 

Thus, to validate the SDAR approximation, we replace "the detailed implementation of the IEEE 802.11 DCF in 
NS-2 MAC layer" with "the SDAR model of contention" keeping all other layers unchanged. The modifications are 
summarized in Section|6] In Section|7]we compare the simulation results obtained from (a) the unmodified NS-2 with 
(b) the SDAR approximation in NS-2, and show that the simulation results obtained from (a) and (b) match extremely 
well. This validates the SDAR approximation for infinite or finite buffers, and equal or unequal arrival rates. In 
Section 15731 we obtain numerical predictions for the performance measures such as collision probability, throughput 
and mean packet delay for the case of equal arrival rates, equal PHY rates, and finite and equal buffers. To proceed 
towards that goal, we first propose a state space reduction technique in Section|4] 

4. Reduction of the State Space 

Suppose that the arrival rates into the queues are equal, i.e., let Aj = A, Vi, 1 < i < M. (We still retain the "infinite 
buffer" assumption in this section.) Thus, Vi, 1 < i < M, Vj, j = 0, 1, . . ., we have di(j) — d(j), Si(j) — s(j), and 
Ci(j) = c(j) where d(j), s(j), and c(j) can be obtained by substituting Aj = A in Equations (O-©. 

Definition 4.1 (Exchangeability 01211 '). The random variables X\ , X 2 , ■ ■ ■ , X n are said to be exchangeable if all the n\ 
permutations (X^ , Xi 2 , . . . , XiJ) of (X\, X 2 , . ■ . , X n ), where i 2 , . ■ ■ , i n ) denotes a permutaion of (1, 2, . . . , n), 
have the same n-dimensional distribution. ■ 

Since the nodes are homogeneous and they use equal PHY rates, and since packets arrive into the node queues 
according to equal rate Poisson processes, it follows that, Vt > 1, the random variables Qi(t), Q 2 (t), ■ ■ ■ , Qm(*), are 
exchangeable provided that the random variables Qi (0), Q 2 (0), . . . , Qm (0), are exchangeable. We assume through- 
out this paper that the random variables Qi(0), Q2(0), . . . , Qm(0), are exchangeable. We remark that, with unequal 
arrival rates into the node queues or with unequal PHY rate of transmission, the desired symmetry will be lost and the 
random variables Qi(t),Q 2 (t), (t), will not be exchangeable. 

Since, Vt > 0, the random variables Qi(t), Q 2 (t), . . . , Qi\i(t), are exchangeable, we consider the following 
alternative description of the system. We define the state of the system at a channel slot boundary t by 



denotes the number of nodes, other than the tagged node, that are non-empty at t. Thus, by definition, we have 





(8) 



i=2 



Note that, Vt > 0, we have, < M{t) < M - 1, and 



N(t) = l iQl{t)>0} +M(t). 
The state space has now reduced from N M to N X {0, 1, ... , M - 1}. 



(9) 



4 The M = 1 case can be easily analyzed, for example, by the "generating function approach." 
5 Any node can be taken as the tagged node due to symmetry. 
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4.1. Approximating the Process {X(t)} by a Process with a Few Unknown Parameters 

In Appendix [C] we derive the one-step transition probabilities of the process {A'(i)} where we show that the 
transition probabilities are functions of some unknown quantities which appear because of the following reason. To 
model the change in M(t) over one step, i.e., over one channel slot, one requires the probability that "a departure 
from a non-tagged node leaves the queue empty." Clearly, a departure leaves the queue empty if the following two 
events occur together: 

(El) In the beginning of the success channel slot, the queue contains exactly one packet (which departs at the end of 
the success channel slot), and 

(E2) The queue does not receive any packets in the success channel slot. 

Event (E2) occurs with probability s(0) (see Equation (0 and recall that the arrival rates are equal). However, with 
the state description X(t), the probability that event (El) occurs cannot be determined for the non-tagged queues, since 
X(t) does not keep track of the number of packets in the non-tagged queues. As shown in Appendix|C] the unknown 
quantities are precisely the probabilities q(i, n, t), t > 0, i = 0, 1, 2, . . ., < n < M — 1, given by 

q(i,n,t) := P(Qi(t) = 1 | Qx(t) = i, M(t) = n, L(t + 1) = L succ , D^t + I) = l,Qi(t) > 0) 

where I, 2 < I < M, denotes the index of a non-tagged node. Since the nodes are homogeneous, and the arrival rates 
and the PHY rates are equal across the node indices, the index I does not really matter, and hence, does not appear in 
the notation q(i, n, t). The condition L(t + 1) = L succ indicates that there is a departure in the (t + l) th channel slot. 
The condition Di(t + 1) = 1 indicates that the departure occured from the I th queue. The condition Qi(t) > must 
accompany the condition Di(t + 1) = 1, since otherwise (Qi(t) = 0) (Di(t + 1) = 0), which is a contradiction. 
Clearly, q(i, n, t) is the probability that the non-tagged queue from which a departure occurs in the (t + l) th channel 
slot contains exactly one packet at t, given the state description (Qi(t) = i, M(t) = n). 

The q(i, n, £)'s are not known precisely because the number of packets in the non-tagged queues are not kept in 
the state description X(t). Moreover, the q(i, n, t)'s cannot be obtained from the known quantities, namely, the arrival 
rate A and the state-dependent attempt probabilities /3„'s. All the other probabilities, namely, Pidie,n, Pcoii.n, Psucc.n, 
< n < M, d(j), c(j), s(j), j — 0, 1, 2, . . ., that appear in the transition probabilities of the process {X(t)} (see 
Equations (IC.27t and (lC.37b ) can be obtained from A and the /3 n 's. 

The conditions L(t + 1) = L succ and Di(t + 1) = 1 can be eliminated from the definition of the q(i, n, i)'s by 
the following lemma. 

Lemma 4.1. Vj > 1, VI, 2 < I < M, Vi, i = 0, 1, 2, . . ., Vn, < n < M - 1, the following is true: 

P(Qi(t)=j | Q 1 {t)=i,M{t) = n,L{t + l) = L succ ,D l (t + l) - l,Q,(t) > 0) 
= P(Qi(t)=j | Qi(t) = i,M(t) = n,Qi(t)>0). 

Proof: See AppendixlBl ■ 
Lemma |4~T1 says that, the probability that a non-tagged queue contains exactly j packets (j > 1), given that it is 
non-empty and given the state description, does not depend on whether a departure occurs from that queue. Thus, 
applying LemmagJ] we redefine, Vt > 0, Vi, i = 0, 1, 2, . . ., Vn, < n < M - 1, VI, 2 < I < M, 

q(i,n,t) :=P(Qi(t) = 1 | Q x {t) = i, M{t) = n, Q,(t) > 0). (10) 

We now apply an approximation first introduced in IBOll in the context of ALOHA networks, and later, also applied 
in Q] in the context of 802. 1 1 WLANs. 

Approximation 4.1 (Conditional Independence). Vt > 0, V£, i = 0, 1, 2, . . ., Vn, < n < M - 1, VI, 2 < I < M, 
we impose the following: 

P(Qi{t) = 1 | Qi(t) = i,M(t) = n, Qi{t) > 0) 
= P(Qi(t) = 1 | Qi{t) = i, N(t) = l {i>0} + n, Qi(t) > 0) 

= P(Qt(t) = l\N(t) = l {i>0} + n,Qi(t)>0). (11) 
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Approximation l4~Tl pertains to the second step in Equation ( TTTT i. which amounts to saying that, the probability that 
a non-tagged queue contains exactly one packet, given that it is non-empty and given the number of non-empty nodes 
in the system, is independent of the exact number of packets in the tagged queue. 

We define, Vf > 0, Vn, 1 < n < M, VI, 2 < I < M, 

q(n,t) :=P(Qi(t) = 1 | N(t)=n,Qi(t) > 0). 
Then, Approximation PTTl savs that, 

M) = \ ?(n+l,t) if *>0. (12) 

We approximate the process {X(t), t > 0} by a process {X(t), t > 0} as follows. The process {X(t),t > 0} has 
the same state description as that of {X(t),t > 0}. However, the transition probabilities of the process {X(t), t > 0} 
are obtained from the transition probabilities of the process {X(t),t > 0} by first applying Lemma |4~T1 and then 
applying Approximation |4~T1 i.e., the transition probabilities of the process {X(t),t > 0} now involve the unknowns 
q(n, t), 1 < n < M. We further impose the conditions that, Vn, 1 < n < M, Vt > 0, 

q(n,t) = q(n), 

where the q(n)'s are unknown constants independent of t. Thus, the process {X(t),t > 0} models the event (El) 
through constant time-independent probabilities q(n)'s. 

With the q(n)'s thus defined, the probability that "a departure from a non-tagged queue leaves the queue empty 
given that it is non-empty in the beginning of the success channel slot and that there are n non-empty nodes in the 
system in the beginning of the success channel slot" is given by q(n)s(0) which is the joint probability of the events 
(El) and (E2). Thus, we can regard the process {X(t), t > 0} as a DTMC (embedded at the channel slot boundaries) 
whose transition probabilities are functions of the unknown parameters q(n) 's yet to be determined. 

In the remainder of this paper, a random variable X (resp. a quantity x) defined for the process {X(t),t > 0} will 
have an analogous random variable X (resp. a quantity x) defined for the process {X(t), t > 0}, and vice versa. Also, 
random variables and probabilities without the time argument would correspond to the stationary regime assuming 
that the stationary regime exists. 

4. 2. Transition Probability Matrix of the DTMC { X (t ) } 

We define the transition probabilities of the DTMC {X(t),t > 0} as follows: 

Aj (n, k) := transition probability from the state (0, n) to the state (j, k), j = 0, 1, 2, . . ., < n, k < M — 1. 
Bj(n, k) := transition probability from the state (i, n) to the state (i + j, k), i = 1,2, . . ., j = —1, 0, 1, . . ., 

< n, k < M - 1. 

The transition probabilities of the process { X(t), t > 0} have been derived in |D](see Equations ( ID. Il l and ( 1D.2| |). 
It follows from Equations ( ID.U and iD.2\ that each of the Aj(n, fc)'s and the Bj(n, fc)'s can be separated into two 
parts, a part that contains an unknown parameter q(-) and a part that does not, and we define 

Aj(n,k) := Af\n,k) + q{n)Af\n,k) } 

B 3 (n,k) := B J (0) (n,fc) + ,?(n + l)B J (1) (n,fc), (13) 

where Af\n,k), Af\n,k), Bf\n,k) and Bf\n,k) can be obtained from Equations dD.ll ) and dD.21 >. and are 
given by 

Afin^k) = ( M ~_ n n _1 ) (Pidle, n d(j) (l-d(0)) k - n d(0Y 

+PcollnC(j) (l-cm^ciOf 1 -*- 1 

+Psucc,n s{j) (1 - S (0)) fe -" s(Q) M - k - 1 ) , (14) 
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Af\n, k) = Psucc , n s(j) (1 - s(0)) K -" S (0) il 



M - n 
k — n + 1 



(l-*(0)) 



M-n-1 

k — n 



(15) 



Bf\n,k) 



M-n-1 
k — n 



M — n- 1\ p. 
k — n 



) (pidle,n+l d(j) {l-d{Q)) k - n d{Q) M - k - 1 

+Pcoll,n + l C(j) (1 - c(0)) k - n ciO)^- 1 

+Psucc,n+1 S(j) (1 - .S(0)) fc -" S (0) M - fc -^ 

) (1 - s(0)) k - n siO)^" 1 (s(j + 1) - s(j)) , 

J n+l 



(16) 



B (1) (n,fc) 



n+l 



" Vcc,„+i {l-a{0)) k - n a(Q) M - k - 1 



M -n 
k-n+1 



(l-s(O)) 



M-n-1 
k — n 



(17) 



Let Aj (0) , A,- (1) , Bj (0) and denote the Af x Af matrices with their (n, fc) t/l entries given by Ay\n, k), 



U o(l) 



Aj (n, fc), Sj- (n,k) and flj (n, fe), respectively. Using this matrix notation, Equation dT3l l can be rewritten as 

Aj := Aj® + A^aAjW , Bj := b/°) + A^B/ 1 ', (18) 

where A^a = diag (0, g(l), . . . , q(M — 1)) and A^b = dia<? (g(l), g(2), . . . , <?(Af )) are Af x M diagonal matri- 
ces. Let jr(j, k), j > 0, < k < M — 1, denote the stationary distribution of the DTMC {X(t), t > 0} (assuming 
that the stationary distribution exists). We define, Vj, j > 0, 



*i := (*(j>0).*0'»l).---.7rCj,M-l)), 



and 



7T :— (7To, 7Ti, 7T 2 , . . . j. 

Using the above notation, the balance equations for the ff (j, k)'s can be written as 

where the transition probability matrix P has the following "M/G/l type" l28ll structure 

A A 1 A 2 A 3 

-B_i Bo Bx B2 

£_! Bo Si 

B_i B 



(19) 



(20) 



5. The Finite Buffer Case 

Let us now turn to the question of determining the unknown parameters <j(n)'s. Applying exchangeability, we can 
write 



q(n,t) := P(Qi(t) = 1 | N(t) = n, Qi(t) > 0) 
= P(Qi(t) = 1 I N(t) = n, Qi(t) > 0) 



(2 < Z < Af) 
(exchangeability). 



(21) 
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Equation (|2TT > provides the clue to estimate the g(n)'s from the stationary distribution of the DTMC {X(t), t > 0} 
as follows: ( recall that random variables and probabilities without a time argument indicate the stationary values) 



q(n) 



p(Qi 
P(Qi 



1\N 
■ 1,N = 



i,Qi >0) 
) P(Qi 



l,M 



P(Qi >0,N = n 
7r(l,n-l) 



P(Qi >0,.M 



(22) 



where recall that 7r(j, fe), j = 0, 1, 2, . . ., < k < M - 1, denote the stationary distribution of the DTMC > 
0}, assuming that it exists. This suggests an iterative method of solution as follows. Given the arrival rate A and 
the state-dependent attempt probabilities /3 ra 's, one can begin with some "guess" values for the q(n)'s and obtain the 
transition probabilities. Then, the stationary state probabilities can be computed and new estimates for the q(n)'s 
can be computed from the stationary distribution by applying Equation (l22l . This procedure can be repeated several 
times until the solutions converge within some specified tolerance limit. However, with infinite buffers, the DTMC 
{ X(t), t > 0} would have infinite number of states and the stationary distribution cannot be computed numerically. 
Thus, we will apply the above idea of estimating the q(n)'s only for the case when the nodes have finite and equal 
buffers so that the finite number of balance equations can be solved numerically. 

Let the buffer size of each queue be K packets (one for the packet in service, if any, and K— 1 waiting for service in 
the queue). We denote the finite buffer version of the process {X{t), t > 0} by {X^ (t),t > 0}. (Processes, random 
variables and probabilities pertaining to the finite buffer case will be denoted by adding a superscript W.) Recall 
that we interpreted the process {X(t),t > 0} as a DTMC embedded at the channel slot boundaries whose transition 
probabilities are functions of the unknown parameters g(n)'s. Similarly, we interpret the process {X( K >(t),t > 0} 
as a DTMC embedded at the channel slot boundaries whose transition probabilities are functions of the unknown 
parameters q^ (n)'s. Let ffW {j, k), < j < K, < k < M - 1, denote the stationary distribution of the DTMC 
{X( K >(t), t > 0} (assuming that the stationary distribution exists). 

We define, Vj, < j < K, 



U,M-1)), 



and 



-CRT) ~(K) ~{K) 



'■— (^0 > "1 > "2 

Then the balance equations for the finite buffer case can be written as 

where the transition probability matrix p(*0 has the following M/G/1/ K structure 



p( K ) = 



(23) 



" A 




A 2 • 








Bo 


B ± 


Bk-2 









Bo 


Bks 










B-i 


Bk-4 


J2j=K-3 Bj 











■ B_i 





(24) 



The Aj's and the Bj's in p( K ) are given by Equation ( fTSI l. except that A,^ and A^b need to be replaced 
by A^(k) >a and A^(ic) >B , respectively, where A^k)^ a = diag (0, gW(l), . . . ,cf K \M — 1)) and A^(ic> iB = 

dw 5 (q (K) (l),g w (2),.--,g (K) ( M )) are M x M diagonal matrices, and A,- (0) , B/ ) and are as 

defined for the infinite buffer case. We emphasize that, p( K ) denotes the transition probability matrix for the finite 
buffer case, and that p( K ) is not the K-step transition probability matrix associated with the one-step transition 
probability matrix P for the infinite buffer case. We also emphasize that computing the infinite sums in the last column 
of p( K ) only requires summing up probabilities from Poisson distributions which can be simplified by observing that 

EJlfc d U) = 1 - E)=o d (i)> and so on - 
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5.1. An Iterative Method 



In this subsection we propose an iterative method which can be applied in the "equal arrival rates, and finite 
and equal buffers" case to obtain the gW(n)'s and the n^ K \j, fc)'s. Given the ^ K \j, /c)'s, we can "estimate" the 
g( K '(n)'s as follows: ( recall that random variables and probabilities without a time argument indicate the stationary 
values) 



gW(n) « P(Q[ K) =l\Q{ K) >0,NW=n) 

-1) 



(25) 



For a given A, assuming some values for the q( K \nys, (f K \n) <E [0, 1], 1 < n < M, the finite number of balance 
equations given by Equation ( 1231 can be numerically solved along with the normalization equation 

K Af-1 

^^^(j,fc) = l (26) 

j=0 fc=0 

to obtain unique positive solutions for the tv^ k > (j, k)'s (provided that {X^ K > (t),t > 0} is positive recurrent for the 
values assumed for the gW(m)'s), Equation ( f25l ) then provides estimates for the q( K \n)'s. The balance equations 
given by Equation ( 1231 and the normalization equation, Equation ( f26b . can be solved with the new set of q( K \nYs 
yielding new estimates for the n^ K \j,kys. This iterative procedure should continue until the solution converges 
within some specified tolerance limit. 

Clearly, for A > 0, the finite-state DTMC {X^ K > (t),t > 0} is irreducible if and only if the q( K \n,ys are positive, 
since Ai^- K \t) can decrease with positive probability if and only if the q( K \nys are positive. Thus, for A > 0, if 
the q^ K '(n)'s are positive, then the DTMC {X( K >(t),t > 0} is irreducible, and hence, it is positive recurrent since 
it has finite number of states. This implies that, if the iterations began with positive q^ K '(n)'s, then, Equations d23l 
and (1261 would yield unique positive solutions for the TX^ K \j, k)'s. Equation ( l25l l would then provide new estimates 
for the cf K> ) (n)'s which would also be positive and this procedure would continue. 

Recall that for a given A, the d{j)\, c(j)'s and s(j)'s are known for all j > (see Equations (O-©) and the j3 n 's 
are obtainable from a saturation analysis ifsl Eoll . Thus, for a given A, the matrices Aj^°\ A^-^ 1 ^, j > 0, and Bj^ \ 
Bj^K j > — 1, need to be computed only once. Each iteration gives anew set of q^ K \n)'s, and the matrices A^K) rA 
and A^k) B get updated in each iteration resulting in updated Aj's and the Bj's (see Equation (fTSTl). This, in turn, 
updates p( K ) resulting in a new set of balance equations which yields a new set of q( K > (n)'s. Important performance 
measures can be obtained from the converged values of the (j, k)'s (Section l53l l. 



5.2. Complexity of our Iterative Method 

The SDAR approximation has the following advantages: 

1. It enables computation of the j3 n 's using the saturation analysis of or ll20ll . Thus, in contrast to the model 
in 1H4I1 . the /3 ra 's in our model are independent of the arrival rates and of the average queue occupancies. They 
are computed only once in the beginning of the overall computation by a separate procedure and then used as 
given parameters in the iterative method described in Section l5TI 

2. The effects of the backoff parameters are effectively captured in the pre-computed /3 n 's (which will be evident 
from the accuracy of our model as demonstrated in Section|7]i. Hence, they are not considered when analyzing 
the queueing dynamics. This enables us to eliminate the first dimension, namely, "the backoff stage of the 
tagged node" from the three-dimensional Markov chain of [ 14] and our model requires computations involving 
a two-dimensional Markov chain. 

3. Since the j3 n 's are computed independent of the arrival rate of the Poisson arrival processes, they need not be 
computed for each arrival rate when studying the effect of arrival rate on performance measures. 
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Thus, the SDAR approximation makes our iterative method computationally less expensive than that in 111411 as 
follows. The complexity of the model in [ 14] is O(RKM) where R denotes the retransmit limit, K denotes the buffer 
size and M denotes the number of nodes (see Section III-C of M14I1 V The complexity of obtaining (3 n ,l < n < M, by 
a separate procedure is O(RM) and the complexity of our finite buffer model, given the /3 n 's, is O(KM). Thus, the 
overall complexity of our finite buffer model is 0(KM) + 0(RM) which is less than the complexity O(RKAl) of the 
finite buffer model of II 1411 . If one needs to solve for, say, I different arrival rates, to examine how the protocol behaves 
with the variation of traffic intensity, then the complexity of our model is 0(lKM) + O(RM) since we compute 
the f3 n '& only once and use them for all the I arrival rates whereas the complexity of the model in lfl4Tl for I different 
arrival rates is 0(lRKM). Thus, when studying the effect of arrival rates on the performance measures, our model is 
far superior to that in lfl4Tl . This reduction in complexity is achieved precisely due to the SDAR approximation. 

5.3. Prediction of Performance Measures 

The processes {Q(t), t > 0} and its finite buffer version {Q^ K ' (t),t > 0} are DTMCs embedded at the channel 
slot boundaries. It is easy to see that {(Q(t),T(t)),t > 0} and {{Q {K) (*), T(t)), t > 0} are Markov renewal 
sequences [191. We apply Markov regenerative analysis to predict performance measures such as collision probability 
and throughput. Also, we obtain the mean packet delay by applying the M/G/l/K queue analysis developed in ll 811 
(see pages 201-205 of d). 

Let {N^Xt)^ > 0} denote the finite buffer version of the process {N(t),t > 0}. Let {N^ K \t),t > 0} be 

derived from {X {K) (t),t> 0} by applying iV (K) (i) := Q[ K) (t) + M (K) (i). We denote the stationary distribution 
of the process {A^ K )(i),i > 0} by p (K} {n), < n < M, and that of the process {N^ K \t),t > 0} by p^^n), 
< n < M. The p( K ^ (n)'s provide "approximations" for the (n)'s and can be obtained from the (j, fc)'s as 
follows: 

p(*)(n) := P(iVW =n)=P{Q[ K) =0,M^ =n)+P{Q[ K) >0,M™ = n - l) 

K 

= n( K 1(0,n) + Y,K (K) (3,n-l). (27) 

3=1 

The TT^ K \j, fc)'s, in turn, can be obtained by the iterative method discussed in Section IBTTl Once the p^^nj's are 
approximated by the p( K> > (n)'s, collision probability and throughput can be obtained as follows. 

Collision Probability: Let A^ K \t) and C^ K \t) denote the total number of attempts and collisions, respectively, 
up to time t, where we recall that t is the discrete time index. The (conditional) collision probability 7 is given by 

M M 

£pW(n)E„(C) £p-W(n)E„(C) 
7 := t^c #f = H w H ' (28) 

n=0 n=0 

where E„(A) and E„(C) respectively denote the mean number of attempts and collisions per channel slot given that 
n nodes are non-empty in the beginning of the channel slot. It is easy to see that 

E n (A) = n(3 n , and E n (C) = n/3 n (l - (1 - M^ 1 ) . 

Throughput: Let (t) denote the total number of successful transmissions up to time t. The aggregate system 
throughput 9 in packets/sec is given by 

M M 

c«w E^ K) («)E„(5) ^pWHE„(5) 
9 := li m !S2 « (29) 

t^oo t M M 

Y,P {K \n)E n {L) ^pWWE n (L) 
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where E n (S) denotes the mean number of successful transmissions per channel slot and E„(i) denotes the mean 
duration of a channel slot in seconds given that n nodes are non-empty in the beginning of the channel slot. It is easy 
to see that 

E n (5')=n ( n (l-/3„r- 1 , 

and 

^-n\L) — \Pidle,n& H~ Pcoll,n 

(a + T c ) + 

Psucc,n 

(<r + T s )) = (a+pcoii 

Due to symmetry, the per-node throughput 6 is given by 

M 

Mean Packet Delay: Consider the process {Q (K) (t),t > 0}. Clearly, for A > 0, {Q {K) {t),t > 0} is an irre- 
ducible and aperiodic DTMC with finite state space (see Appendix lAl where we prove irreducibility and aperiodicity 
of {Q(t),t > 0}). Hence, for A > 0, the DTMC {Q^ (t), t > 0} is positive recurrent, stationary and ergodic. For 
any tagged queue, define the following: 

a(j) := fraction of (real) time that the tagged queue contains j packets, < j < K. 

P (j) '■= probability that a departure from the tagged queue leaves j packets behind, < j < K — 1. 
(Note that the number of packets left by a departure cannot be more than K — 1) 

p( a \j) := probability that a packet accepted into the tagged queue finds j packets, < j < K — 1. 

(Note that the number of packets as seen by an accepted packet cannot be more than K ~ 1) 

In Appendix[E] we obtain the probabilities p^ d \j), < j < K — 2, as follows: 

j+l M-l , s 

p {d) U) = - n m-! : ■ (3Q) 

,;=1 n =a ^ 

where 7r^^(j,n), < i < K, < n < M — 1, denotes the stationary probability that the tagged queue contains i 
packets and n of the non-tagged queues are non-empty. The s(j — i + l)'s in Equation (f30b can be obtained using 
Equation (0 and the (i, n)'s can be approximated by the (i, n)'s. The (i, n)'s, in turn, can be obtained 
by the iterative method discussed in Section l5TI The probability (K — 1) is given by 

K-2 

p^(K-l) = l-J2p {d) U)- (3D 

Since arrivals and departures occur one at a time, a level-crossing analysis gives, Vj, < j < K — 1, 

pW(j)=p(«0(j). (32) 

Note that the probability that an arrival is blocked is given by a(K). The mean rate at which packets are accepted 
into the queue must be equal to the mean rate at which packets depart from the queue. Hence, a(K) can be obtained 
by solving 

X(l-a(K)) = e = ^, (33) 



where 9 can be computed applying Equation ( 1291 ). According to [18], we have, Vj, < j < K — 1, 



"(i) =p [a \j){l-a{K)). (34) 
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The mean queue length Q is then given by 



K 



(35) 



3=0 



from which the mean sojourn time or the mean packet delay W can be obtained as 




(36) 



6. Applying the SDAR Approximation in the NS-2 Simulator 

Wireless network simulators invariably employ simple propagation models at the PHY layer to keep the simu- 
lations reasonably fast. In this section, we describe a model-based simulation technique at the MAC layer which is 
based on the SDAR approximation. We apply the SDAR approximation in NS-2 as follows. 

We modify the original implementation of NS-2 MAC layer to keep track of the number of non-empty nodes. Ar- 
rivals that occur during activity periods (i.e., during successful transmissions or collisions) are not taken into account 
to update the number of non-empty nodes until the activity finishes. Hence, the number of non-empty nodes does not 
change during the activity periods. Whenever an activity finishes or an arrival occurs during a channel idle period, the 
number of non-empty nodes is updated. Whenever the number of non-empty nodes is updated, all previously sched- 
uled transmissions (if any) are canceled and random backoffs are sampled for each non-empty node using independent 
geometric random variables each having a mean 1 / (3 n where n denotes the current value of the number of non-empty 
nodes. Note that the geometric backoff durations with mean l/j3 n are equivalent to Bernoulli attempt processes with 
probabilities /3„. The /3 n 's are pre-computed using the saturation analysis in f20ll and stored in a look-up table. 

From the sampled backoffs, it is easy to determine which node(s) sample the minimum backoff. If only one node 
samples the minimum backoff, the next event is a successful transmission. If two or more nodes sample the same 
minimum backoff, the next event is a collision. The appropriate event is then scheduled. If arrivals occur to empty 
queues before the "beginning of the scheduled event" epoch to increase the number of non-empty nodes, then the 
scheduled event is canceled. Else, the scheduler clock is moved to the "end of the scheduled event" epoch. In case 
the next event is a successful transmission, the DATA frame is handed over to the destination's MAC layer which then 
generates the corresponding ACK frame. 

The above modifications have enabled us to achieve speed-ups up to 5.4. These speed-ups have been achieved with 
respect to the MAC layer operations and are summarized in Tables Q]and[2]for two different machines. The observed 
speed-ups are obtained due to the following reasons. In NS-2, one transmission event per non-empty node remains 
pending in the "scheduler queue" and n backoff timers are kept running for n non-empty nodes. When a timer expires 
resulting in a transmission, all the other n — 1 timers remain "paused" until the transmission finishes. When the timers 
resume, the remaining n — 1 transmission events have to be rescheduled. Similar pausing and rescheduling occurs in 
case of collisions when multiple timers expire simultaneously. Note that rescheduling of events requires searching the 
appropriate events from the scheduler's event queue which consumes the largest fraction of MAC simulation time. In 
our modifications, due to the memoryless property of the SDAR attempt model, we do not have to keep the backoffs 
sampled by the nodes. Moreover, at any point of time, only one event remains pending at the MAC layer which is the 
next event to occur on the common channel. This single pending event is interrupted and is rescheduled only if an 
arrival increases the number of non-empty nodes. Hence, the speed-up increases with the number of nodes M. The 
speed-up also increases with arrival rate A since the average number of non-empty nodes increases with increase in 
A. However, above a certain A, the rate of cancellation of already scheduled events dominates and speed-up actually 
decreases with A. The speed-up w.r.t. MAC layer operations remains constant in saturation conditions even though 
the simulation time consumed in higher layer operations increases with A due to increasing number of arrivals. 

These observations are supported by the data in Tables Q] and [2] where the last row for each M corresponds 
to saturation. The increase in speed-up with M is particularly desirable since NS-2 is found to become worse with 
increase in M than A. Also, note in Tables[T]and|2]that the speed-ups are more for the faster machine, i.e., MACHINE- 
I. This indicates that the speed-ups are not due to the incapability of the machines. 
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Table 1: Speed-up for MACHINE-I (Pentium dual core, 2.80 GHz, 1024 KB cache) 



M 


A 

(pkts/sec) 


MAC time 
NS-2 (sec) 


MAC time 
SDAR (sec) 


Speed-up 


50 


5 


45.045 


9.405 


4.79 


50 


10 


103.43 


19.16 


5.4 


50 


15 


158.44 


39.40 


4.02 


30 


5 


10.8 


3.25 


3.08 


30 


10 


22.88 


6.89 


3.32 


30 


15 


32.49 


10.24 


3.17 


30 


20 


60.49 


22.84 


2.65 


10 


10 


2.233 


0.943 


2.37 


10 


20 


4.314 


1.864 


2.31 


10 


30 


7.256 


3.114 


2.33 


10 


40 


8.351 


3.641 


2.30 


10 


50 


10.777 


4.767 


2.26 


10 


60 


13.269 


5.789 


2.30 


10 


70 


19.922 


10.312 


1.93 



Table 2: Speed-up for MACHINE-II (Pentium, 1500 MHz, 256 KB cache) 



M 


A 

(pkts/sec) 


MAC time 
NS-2 (sec) 


MAC time 
SDAR (sec) 


Speed-up 


50 


5 


102.06 


31.99 


3.19 


50 


10 


268.23 


66.86 


4.0 


50 


15 


369.13 


97.77 


3.78 


30 


5 


23.9 


9.59 


2.49 


30 


10 


48.16 


19.06 


2.53 


30 


15 


75.93 


29.79 


2.55 


30 


20 


136.82 


54.11 


2.53 


10 


10 


3.66 


2.13 


1.72 


10 


20 


7.31 


4.25 


1.72 


10 


30 


11.38 


6.62 


1.72 


10 


40 


15.52 


9.00 


1.72 


10 


50 


20.25 


11.28 


1.795 


10 


60 


25.87 


13.99 


1.85 


10 


70 


31.07 


20.02 


1.55 



7. Numerical and Simulation Results 

In this section we provide numerical and simulation results to validate our coupled queue model and our iterative 
method. The values of the backoff parameters were taken as per the 802.11b standard. We took Basic Rate = 2 
Mbps, Data Rate =11 Mbps, Packet Payload Size = 1000 Bytes. We provide results for the Basic Access mode. 
The iterative method was implemented using MATLAB. All other protocol parameters are set as per the standard 
setting in 802.11b. Figures [5] [6] andQcompare the collision probability 7, the per-node throughput 9, and the mean 
packet delay W with M = 10 nodes, infinite buffer size, and equal arrival rates. The buffer sizes were set to very 
large values to simulate infinite buffers. It can be seen that the simulation results obtained from the unmodified NS- 
2 and the SDAR approximation in NS-2 match extremely well. The mismatch in the collision probabilities which 
leads to visible mismatch in the mean packet delays near saturation is mainly due to the over-estimation of collision 
probability by the saturation analysis H, |2(J, which clearly appears as a w 5% mismatch of collision probability 
beyond the saturation threshold (which is about A = 62.5 packets/sec for M = 10 nodes, (see Figure|5]l). 

Notice, however, in Figure|6]that the throughputs obtained from NS-2 simulations and from the SDAR approxima- 
tion in NS-2 are indistinguishable. Let <d sa t denote the aggregate system throughput under saturated conditions. For 
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Figure 5: Comparison of collision probability 7 with M = 10 
nodes, infinite buffer size, and equal arrival rates. 



Figure 6: Comparison of throughput per node 9 with M = 10 
nodes, infinite buffer size, and equal arrival rates. 
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Figure 7: Comparison of mean packet delay W with M = 10 
nodes, infinite buffer size, and equal arrival rates. 



Figure 8: Comparison of collision probability 7 with M = 10 
nodes, finite buffer size K = 5, and equal arrival rates. 
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Figure 9: Comparison of throughput per node 6 with M = 10 
nodes, finite buffer size K = 5, and equal arrival rates. 
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Figure 10: Comparison of mean packet delay W with M = 10 
nodes, finite buffer size A" = 5, and equal arrival rates. 



the default backoff parameters as prescribed in the 802.1 lb standard, O sat is not sensitive to the change in . Hence, a 
less than 5% difference in 7 does not lead to noticeable change in 8 sa t , and Q sa t in the original NS-2 is almost equal 
to Q sa t in NS-2 with SDAR. Since the number of nodes, M, remains constant, 9 sat := Q sa t/M, in the original NS-2 
is also equal to 9 sat in NS-2 with SDAR. This explains the good match of the flat portion of the curves in saturation. 
Since there are no packet drops due to buffer overflows (infinite buffer) and there are no packet discards, one can 
expect that, for A < 9 sa t, the rate of departure is equal to the rate of arrival, i.e., 9 = A, since the queues are stable. 
This explains the "y = x" behavior below saturation load. 

Figures [8] |9] and [TO] compare the collision probability 7, the throughput per node 9 and the mean packet delay W 
with M — 10 nodes, a buffer size of K — 5 packets per queue, and equal arrival rates. It can be seen that the SDAR 
approximation in NS-2 and our iterative method of analysis both match extremely well with the unmodified NS-2. 
Furthermore, the results from the SDAR approximation in NS-2 and that from our iterative method of analysis are 
extremely close. This validates the conditional independence approximation, Approximation |4~T1 (at least from the 
point of view of predicting the performance measures), and also validates our iterative method. 

Notice the peak of throughput in Figure [9] which occurs primarily because of the finite buffer size, and thus, not 
visible in the infinite buffer case. Note that, in the finite buffer case, the system cannot saturate for finite arrival rates. 
This can be understood by considering the finite state DTMC {Q^ K '(t), t > 0}, which is positive recurrent for all 
finite arrival rates, < A < 00. Hence, every queue must be empty for a positive fraction of time. However, as the 
arrival rate A increases, the fraction of time for which a node remains empty decreases. For very large arrival rates 
one can expect the throughput to become close to the saturation throughput. 

The peak throughput in Figure|9]is slightly higher than the saturation throughput because of the following reason. 
Consider an arrival rate A higher than the per-node saturation throughput. Had the buffer sizes of the nodes been 
infinite, each node would have been saturated, the number of contending nodes would have been M, and the aggregate 
system throughput would have been <d S at,M ■ However, due to finite buffers, the effective number of contending nodes 
n e ff(X) corresponding to A, is smaller than M. Since n e ff(X) < M, we have O sa t,n e// (A) > Qsat,M, and the 
aggregate system throughput is larger than the saturation throughput. With increasing A, n e j f ( A) would approach M 
and the aggregate system throughput would approach <d S at.M- 

We now turn to some simulation results with unequal arrival rates. Figures [TT1[l"4l compare the collision probability 
7, the throughput, the mean service timef^ and the mean packet delay W with M = 2 nodes, infinite buffer size, and 
unequal arrival rates with ratio Ai : A2 = 1 : 2. Figures [151(181 compare the collision probability 7, the throughput, 
the mean service time, and the mean packet delay W with M = 3 nodes, infinite buffer size, and unequal arrival 
rates with ratio Ai : A2 : A3 = 1 : 2 : 3. Notice that in Figures [711181 we plot the results against the arrival rate 
Ai into Node-1. The corresponding arrival rate into Node-2 (resp. Node-3) is given by A2 = 2Ai (resp. A3 = 3Ai). 
It can be seen in Figures [TT|[T8l that results from the SDAR approximation in NS-2 match extremely well with the 
unmodified NS-2. This validates the SDAR approximation also for the case when the arrival rates are unequal and 
differ significantly from each other. 

Notice in Figure [TT] that 7 for Node-1 is larger than that for Node-2 even though Ai < A2 — 2X±. This can 
be explained as follows, recall that 7 is the conditional collision probability, i.e., 7 is the probability with which 
an attempted transmission collides. The probability that "Node-2 is non-empty given that Node-1 is non-empty" is 
larger than the probability that "Node-1 is non-empty given that Node-2 is non-empty" since A2 > Ai. Hence, the 
probability that Node-l's attempts collide is larger than the probability that Node-2's attempts collide. The fact that 7 
for Node-1 is larger than that for Node-2 explains why the mean service time of Node-l's packets is larger than the 
mean service time of Node-2's packets (see Figure [T3l>; a larger 7 implies that a larger number of attempts is need for 
successfully transmitting any packet. However, since A2 > Ai, the "queueing delay" for Node-2 is larger than that for 
Node-1. This explains why the mean packet delay W, which includes the queueing delay, is larger for Node-2 than 
that for Node-1 (see Figure[l4|i. Similar explanations hold for the three-node case as well. 



6 Mean service time is the mean time a packet spends at the Head Of Line (HOL) position until it is dequeued. A packet is dequeued either when 
an ACK is received which implies a successful transmission or when a timeout occurs and the the retransmit limit has been reached. 
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Figure 11: Comparison of collision probability 7 with M = 2 nodes, 
infinite buffer size and unequal arrival rates. 



Figure 12: Comparison of throughput per node 6 with M 
infinite buffer size and unequal arrival rates. 
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Figure 13: Comparison of mean service time with M = 2 nodes, infi- Figure 14: Comparison of mean packet delay W with M = 2 nodes, 
nite buffer size and unequal arrival rates. infinite buffer size and unequal arrival rates. 



8. Conclusion 

In this paper, we applied the SDAR approximation to model the attempt processes of the nodes in a single cell 
WLAN. Applying the SDAR approximation and considering Poisson arrivals, we developed a Markov model and 
reduced the problem of modeling a single cell WLAN with non-saturated nodes to the problem of analyzing a coupled 
queue system. We provided a sufficient condition under which the joint queue length Markov chain is positive recur- 
rent. Assuming equal arrival rates, we proposed a technique to reduce the state space of the coupled queue system. 
For the case of equal arrival rates, and finite and equal buffers, we proposed an iterative method to obtain the station- 
ary distribution of the reduced state process and obtained accurate numerical predictions for important performance 
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Figure 15: Comparison of collision probability 7 with M = 3 nodes, 
infinite buffer size and unequal arrival rates. 



Figure 16: Comparison of throughput per node 6 with M = 3 nodes, 
infinite buffer size and unequal arrival rates. 
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Figure 17: Comparison of mean service time with M = 3 nodes, infi- Figure 18: Comparison of mean packet delay W with M = 3 nodes, 
nite buffer size and unequal arrival rates. infinite buffer size and unequal arrival rates. 



measures. Our iterative method was shown to be better in terms of computational complexity than an earlier work 
with same level of accuracy. The reduction in complexity was achieved precisely due to the SDAR approximation. 
We demonstrated via simulations that the SDAR model of contention provides an accurate model for the CSMA/CA 
protocol (for equal or unequal arrival rates and infinite or finite buffers). In addition, when applied to the NS-2 sim- 
ulator it improves the simulation speed without affecting the accuracy of results. However, a theoretical justification 
for the SDAR approximation is lacking at this point. 
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Appendix 



A. Proof of Theorem [37T1 



We apply Foster's criterion (see page 29, chapter 2 of [ 1 1]) to prove Theorem l3.ll 



Theorem A.l (Foster's Criterion). A time homogeneous irreducible aperiodic Markov chain {X (t), t > 0} with a 
countable state space S is positive recurrent if and only if there exists a non-negative function f(x), x £ S, a number 
e > 0, and a finite set A C S, such that the following conditions hold: 

(1) E(f(X(t + 1)) - f(X(t)) | X(t) = x) < -e, Vx E S \ A, 

(2) E(/(X(t+l)) | X(t) =x) <oo,Vx G A ■ 

Clearly, the DTMC {Q(t),t > 0} is time homogeneous and has a countable state space S = N M . Irreducibility 
and aperiodicity can be easily shown as follows. If, Vi, 1 < i < M, A; > 0, then any state k = (ki, k 2 , . . . , km) S 
N M can be reached from the state = (0, 0, . . . , 0) S N M in one step by ki, 1 < i < M, arrivals to the i th queue. 
Also, the state (0, 0, ... , 0) can be reached from any state (fci, ^2, . . . , fcju) £ N M in YlfLi ki steps by 53i=i ^» 
consecutive successes and no arrivals such that the fc/s do not increase in between. Hence, the DTMC {Q(t),t > 0} 
is irreducible. Since {Q(t),t > 0} is irreducible and there exists a self loop, e.g., from the state (0, 0, ... , 0) to itself, 
the DTMC {Q(t), t > 0} is aperiodic as well. 

We define the finite set A and the non-negative function /(•) as follows: 

M 

A := {0} = {(0, 0, . . . , 0)} ; f(h, k 2 ,..., k M ) := £ h. 

i=l 

Notice that A := {0} = {(0, 0, . . . , 0)} consists only of the "system-empty" state, and for any state fc = (fci, fca, . . . , /cm), 
the function f(ki, k 2 , . . . , fcjv/) := J2i=i ^« gives the total number (of packets) in the system. 

Since the transition probabilities of the DTMC {Q(t),t > 0} in any state depend on the number of non-empty 
nodes in that state, we partition the state space S = N M as 



s = |J s n , 



ri=0 



where 



M 

k = (k u k 2 ,...,k M )eN M : 5> {fc! >o} 



n 



For all k S S n , we have 

E(,f(Q(t + l))-f(Q(t)) | Q(t)=k) 
= E(/(Q(i + 1)) - f(Q(t)) | Q(t) = k, N(t) = 
= P(L(t + 1) = L ulle | N(t) = n) E(f(Q(t + 1)) - f(Q(t)) \ Q(t) = k, N(t) = n, L(t + 1) = L ldle ) 

N{t) = n) E(f(Q(t + 1)) - f(Q(t)) | Q(t) = k, N(t) = n, L(t + 1) = L coll ) 
N(t) = n) E(/(Q(i + 1)) - f{Q{t)) | Q(t) = k, N(t) = n, L(t + 1) = L succ 

- Qi(t)) | Q(t) = k, N{t) = n, L(t + 1) = L ldl )j 
1) - Qi(fif) | Q(t) = k, N{t) = n, L(t + 1) = L col )j 
1) - Qi{t)) | Q(t) = k, N(t) = n, L(t + 1) = L succ ) 
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Vidle.n 

E[5^Ai(t + l) \L(t + l) = L 

idle^J H~~ Pcoll,n 
"I - Psucc,n StlCC I 

Pidle,n iJ2 E (M t + l ) \L(t + l) = L 

idle^J + Pcoll,n coZi / I 



T" Psucc^n succ) I Psucc,n 

\i=l ) 

(M \ / M \ / M \ 

2J AjCT J + p co H, rl I ^ Aj (CT + T c ) J + p SU cc,n I ^ A, (CT + T a ) J - Psucc,n 

( M \ 

= I ^ ] Aj J ((J + Pcoll,nT c + Psucc,nTs) ~ Psucc,n, (A.l) 
\i=l / 

where we have used the following facts: (i) (Q(t) = k 6 5„) => (iV(t) = n), (ii) given 7V(i) = n, the type of (t+l) th 
channel slot does not depend on Q(t), and (iii) given the type of channel slot, the quantity X)f=i {Qi(t + 1) — Qi(t)), 
which represents "the increase in the total number (of packets) in the system in the (t + l) th channel slot," depends 
neither on Q(t) nor on N(t). In particular, (a) if the (t + l) th channel slot is an idle or a collision channel slot, then 
52i=i (Qi(t + 1) — Qi(t)) is equal to Yli=i Ai(t + 1) which represents the total number of arrivals in the (t + l)*' 1 
channel slot, and (b) if the (t + l) th channel slot is a success channel slot, then YliLi (Qi(t + 1) — Qi{t)) is equal to 
{^2iLi A%(t + l)j — 1 since there is exactly one departure in the it + l) th channel slot. 

Noticing that the finite set A = So, and that Pidie.o = l,Pcoll,a = Psuccfl = 0, condition (2) of Theorem lA.il 
is satisfied for the DTMC {Q(t),t > 0} if fe^i A;) cr < oo, i.e., if, Vi, 1 < i < M, X t < oo. Condition (1) of 
Theorem lATTl is satisfied for the DTMC {Q(t), t > 0} if, there exists e > such that, Vn, 1 < n < M, we have 



' M \ 
\i=l } 



' Pcoll^nTc ~t~ Psucc^n^s) Psucc,n — 



Since e can be made arbitrarily small, condition (1) of Theorem lA.ll is satisfied for the DTMC {Q(t),t > 0} if, 
Vn, 1 < n < M, we have 

/ M \ 

Psucc : n 



5> 

vi=l 



< 



((J -j- Pcoll,n-^c Psucc.ri^s) 



P succ n 

We define L sat ^ n := {cr + p coU . n T c + p SUCCyn T s ) and Q sa t,n ■= * CC ' n ■ Note that, L sat . n and 6 sat ,„ represent 

the mean channel slot duration in seconds and the aggregate throughput in packets/sec, respectively, in a single cell 
consisting of n homogeneous and saturated nodes B20I1 . Condition (1) of Theorem lA.ll is satisfied for the DTMC 

{Q(t),t > 0} if, Vn, 1 < n < M, we have 



/ M \ 
\i=l / 



Since the condition (2) of Theorem |AT| holds at all finite arrival rates, the DTMC {Q(i), t > 0} is positive recurrent 
if 



/ m \ 
/ Aj < min Q s 

\^ I l<n<M 



and Theorem l3.1l is proved. 
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B. Proof of LemmaED 

Vj > 1, V7, 2 < I < M, Mi, i = 0, 1, 2, . . ., Vn, < n < M - 1, we have 

P(L(t + 1) = L succ , D t (t + 1) = 1 | Qi(t) = i,-M(t) = n, Q t (t) > 0, Q,(t) = j) 
= P(L(t + l)=L succ | Q 1 {t)=i,M{t)=n,Q l (t) >0) 

■ P(A(i + l) = 1 | Qi(t) = = n,Qi(t) > 0,L(t + l) = L succ ) 

= P(L(t + l) = L succ ,Di(t+l) = l\Q 1 (t) = i,M(t)=n,Qi(t)>6), (B.l) 

since (i) given that Q\{t) — i and Ai(t) — n, we know that N(t) = 1{;>o} + n, (ii) given N(t), the probability 
that L(t + 1) = L succ , does not depend on the exact value j of Qi(t) and is equal to p SU cc,N(t), an d (iii) given N(t), 
Qi(t) > and L(t + 1) = L succ , the probability that Di(t + 1) = 1, does not depend on the exact value j of Qi(t) 
and is given by (see Equation (|7|i)- Notice that, the condition Qi(t) > cannot be eliminated from the second 
factor in the first step (on the right hand side of Equation ( IB.ll i). since Qi(t) = would imply Di(t + 1) = 0. In fact, 
given N(t), the probability that L(t + 1) — L succ does not depend on anything else, and the condition Qi(t) > 
could have been eliminated from the first factor in the first step. However, we keep the condition Qi(t) > since it is 
required in the second factor in the first step, and hence, is required to obtain the final expression. 
Define now the events 

E A := {L(t + I) = L succ , D t {t + I) = 1}, 
E B ■= {Q 1 {t)=i,M{t)=n,Q l {t) > 0}, 
E c :={Q l (t)=j},j>l. 
Equation (IB.U says that, P(E A | E B ,E C ) = P(E A | E B )- Lemma HTT1 follows from 

(P(E A | E B ,E C ) = P(E A | E B )) «■ (P{E C \ E Bl E A ) = P(E C \ E B )). ■ 

C. Transition Probabilities of the Process {X{t)} 

In this appendix we derive the one-step transition probabilities of the process {A'(i)} from a generic state (i, n) 
to a generic state (j, k), i,j = 0, 1, 2, . . ., < n, k < M — 1. 

C.l. When the Tagged Queue is Empty in the Initial State 

Consider the transition probability P{Q\(t + 1) = j, Ai(t+ 1) = k I Qi{t) = Q,A4(t) = nj , which corresponds 
to the case when the tagged node is empty at the channel slot boundary t and can be expanded according to whether 
the (t + l) th channel slot is an idle, a collision and a success channel slot as follows: 

P(Q 1 (t + l)=j,M(t+l) = k | Qtfjt) =0,M(t) = n) 

= P(Q 1 (t + l)=j,M(t+l) = k,L{t + l) = L idU | Q 1 (t)=0,M(t)=n) 
+ P{Qi(t + 1) = j, M(t + 1) = k, L(t + 1) = L coU | Qt(t) = 0,M(t) = n) 
+ P(Qi{t + l) = j,M(t + l) = k,L(t + l) = L succ | Qi(t) = 0,M(t) = n). (C.l) 

C.l.l. The Case of Idle Channel Slot 

The first term on the right hand side of Equation ( 1C.U corresponds to the case when the tagged queue is empty at 
the channel slot boundary t and the (t + l) th channel slot is an idle channel slot, and can be expanded as follows: 

P(Qi(t + l)=j,M(t + l) = k,L(t + l) = L id i e | Qi(t) =0,M(t) = n) 
= P(L(t + l)=L ld i e | Qi(t) =0,M(t) =n) 

■ P(Qi(t + l)=j,M(t + l) = k\ Q 1 {t) = 0,M{t) = n,L(t + l) = L idle ) 
= P(L{t + l)=L ldle | Q 1 (t)=0,M(t)=n) 

■ P(A 1 (t + l) =j,M{t + l) = k | Q 1 {t) = 0,M(t) = n,L(t + l) = L idle ) 
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= P(L(t + l)=L ldle \N(t) = n)- P(M(t + l)=k\Qi{t)=0,M(t)=n,L(t+l) = L idle ) 

■P(A 1 (t + l)=j | L(t + 1) = L ldle ) 
= Pid i e , n -d{j)-P(M(t + l) = k\Q 1 (t) = 0,M(t) = n,L{t + l) = L idle ), (C.2) 

where we have applied Equation © to write P[L{t + 1) = Lidie \ N{t) = nj = Pidie,n- We have also used 
the following facts: (a) (L(t + 1) = L ldle ) => (D m (t + 1) = 0, 1 < m < M), (b) (Qx(t) = 0, D^t + 1) = 0) => 
(Qi(t + 1) = A x (t + 1)) which implies that Q x {t + 1) = j, if and only if A x (t + 1) = j, (c) given that L(t + 1) = 
I/idie, the probability that A x (t + 1) = j, does not depend on anything else and is equal to d(j). (See Equation (fj) 
and recall that the arrival rates are equal). 

The conditional probability P(M{t + 1) = k \ Qi(t) = 0,M(t) = n, L(t + 1) = L ldle ) in Equation (IC.2t can 
be expanded as follows: 



P(M(t + l) = k | Qi(t) = 


0,M(t) 


= n, L(t 


+ 1) — L id i e ) 


P(M(t) =n,M(t + 


- 1) = fc | 


Qi(t) = 


0, L(i + 1) = L id ; e ) 


P(M(t) = n 


| Qi(t) - 


= 0,L(t- 


I- 1) = iid/ e ) 



(C3) 

Let Im '■= {2, . . . , M} denote the set of non-tagged queues. Then, the numerator on the right hand side of 
Equation (IC.3b can be expanded as follows: 

P(M(t) = n,M(t + l) = k\ Qx{t) = 0, L(t + 1) = L idle ) 

J2 p (Qm(t) > 0,Vm e /,Q r (t) =0,VreI M \I,M(t + l) = k\Qi(t) = 0,L(t + l) = L, d / e ) 

/C/m : |/|=n 

^ P{Q m {t) > 0,Vm e /,Q r (t) = 0,Vr e 7 M \ / | QiO) = 0,i(t + l) = L ldle ) 

/C/m : |/|=n 

• P(M(t + 1) = fe I Q m (t) > 0, Vm € /, Q P (t) = 0, Vr e J M \ /, Q x (t) = 0,L(t + 1) = L idie ). 

(C4) 



The second conditional probability inside the summation in the second step of Equation ( IC.4l i can be simplified as 
follows: (Remember that / C Im and |/| = n.) 

P(M(t + 1) = k | Q m (t) > 0, Vm e /, Q P (t) = 0, Vr e / M \ /, Qi(t) = 0, L(t + 1) = L ld/e ) 

= P(M(t + 1) = fc | Q m (i) > 0, Vm e /, Q r (t) = 0, Vr 6 / M \ /, Q^t) = 0, L(f + 1) = L ldle , 

D u (t+l) = Q,\fueI M ) 

= p (\^2 1 {(Q m (t)+A m (t+i))>o}+ ^2 l{A r (t+i)>o} J = k | Q m (t) > 0,Vm e I,L(t + 1) = L id ie) 

\m£l reI M \I J 

= P{ ^ !{A r (t+i)>o} = - n | + 1) = L ldie ) 



M-n-1 
k — n 



{l-d{Q)) k - n d(0) M - k - 1 . (C.5) 



In the first step we have used the fact that (L(t + 1) = Li d i e ) ==> (D u (t + 1) = 0, Vu 6 Im)- In the second step 
we first applied Equations ([8]) and (Q~|) to write 

M M 

M(t+1) = ^ l{Q u (t+l)>0} = ^2 1 {(Qu(t)-D u (t+l)+A u (t+l))>0}, 
u=2 u=2 
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and then substituted D u (t + 1) = 0, 2 < u < M, and Q r (t) = 0, Vr e Im \ I- In the third step we have used the 
fact that Q m {t) > 0, Vm e /, implies that Ylmei ~*-{(Qm(t)+A m (t+i))>o} — \I\ = n - I n tne l ast ste P we observe that, 
SreiM V l{A-(*+i)>o} — — 7i, if and only if arrivals occur to exactly k — n of the \Im \ I\ = M — n—1 non-tagged 
queues that are empty at t, and, given that Lit + 1) = -L^e, the probability that arrivals occur (resp. do not occur) to 
a queue is given by 1 — d(0) (resp. d(0)). 

We recall that by definition, [ ] > 0, only if a > b > 0, and [ ] =0, otherwise. 

. . \ b J . . \ b J 

Substituting Equation dC.5b in Equation dC.4l ). we obtain 

P(M(t) = n,M{t+l) = k | Qi(t) = 0,L(t+l)=L ilfle ) 

P(Q m (t)>0,VmeI,Q r (t)=0,VreI M \l\Qi(t) = 0,L(t+l)=L idle ) 



ICIm ■■ \I\=n 



M — n — 1 
k — n 



{l-d{Q)) k - n d{0) 



M-k-l 



= P(M(t)=n\Q 1 (t) = 0,L(t+l) = L ldle )- ^ ^ ^(1 - d(0)) k ~ n d(p) M - k -\ (C.6) 

where, in the second step we first take the common factor ( M fc "_!7 1 ) (1 - d(0)) k - ,l d{0) M - k - 1 outside the summation, 
and then observe that the summation is equal to the probability P(A4(t) = n Qi(t) = 0, L(t + 1) = Lidie). 
Substituting Equation ( IC.6I ) in Equation ( IC.3I ). we obtain 



P(M(t + l) = k\ Q 1 (t)=0,M(t)=n,L(t+l) = L idle ) = f M fe r ) (1 - ^O))*""^)^- 1 , (C.7) 



where we observe that the probability P(A / ((i) = n Qi(i) = 0, L(t + 1) = Lidie) gets canceled from both the 
numerator and the denominator. 

Substituting Equation dC.7b in Equation dC.2b we obtain 



Pidie,n ■ d(j) ■ ( M , H V ) (1 - d(Q)) k - n d(0) M - k -\ (C.8) 



P(Q 1 (t + l)=j,M(t + l) = k,L(t + l)=L ldle | Q 1 (t)=0,M(t)=n) 

[ — n - 
k — n 

C.1.2. The Case of Collision Channel Slot 

The second term on the right hand side of Equation ( IC.U corresponds to the case when the tagged queue is empty 
at the channel slot boundary t and the (i + l)*' 1 channel slot is a collision channel slot. By similar arguments, we 
obtain 



Pco«,n-c(i)-( M , n l \l-c{0)) k - n c{0) M - k -\ (C.9) 



P(Q 1 (t + l)=j,M{t + l) = k,L(t + l) = L coll | Q 1 (t) = 0,M(t) = n) 

' — n - 
k — n 

C.1.3. The Case of Success Channel Slot 

The third term on the right hand side of Equation dC.U corresponds to the case when the tagged queue is empty at 
the channel slot boundary t and the (t + l) th channel slot is a success channel slot. For this case, we have 

P(Qi{t + 1) - j,M{t + 1) = k, L(t + 1) = L succ | Qi(t) = 0,M(t) - n) 

= P(L(t + l) = L succ | N(t) = n) ■ P(M(t + 1) = fc | Q 1 (t)=0,M(t) = n,L(t+l)=L succ ) 

■P(A 1 (t + l)=j \L{t + l) = L succ ) 
= Psucc,n ■ s(j) ■ P(M(t + 1) = k | Q x {t) = 0,M{t) = n, L(t + 1) = L succ ) . (CIO) 
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Note that, since the tagged queue is empty at the channel slot boundary t, the departure in the success channel slot 
must occur from a non-tagged queue. Thus, we have Q\(t + 1) = j, if and only if A\ (t + 1) = j. Furthermore, given 
that L(t + 1) = L succ , the probability that A\ (i + 1) = j, does not depend on anything else and is equal to s(j). (See 
Equation (01 and recall that the arrival rates are equal). 

The non-tagged queue from which the departure occurs is empty at the channel slot boundary t + 1 if it has exactly 
1 packet at the channel slot boundary t and it does not receive any packets in the (t + \ ) th channel slot. Thus, we must 
condition on the number of packets at the channel slot boundary t in the non-tagged queue from which the departure 
occurs at the channel slot boundary t + 1. To that end, we first expand the conditional probability in the second step 
in Equation ( IC. 10b as follows: 



P(M(t + l)=k | Q 1 (t) = 1 M(t)=n,L(t + l) = L succ ) 

P(M(t)=n,M(t + l) = k | Qi(t) =0,£(i + l) = L s 
P(M(t) = n I Q^t) = 0, L(t + 1) = L succ ) 



(C.ll) 



As before, let 1m '■= {2, . . . , M} denote the set of non-tagged queues. Then, the numerator on the right hand side 
of Equation (IC. 1 U can be expanded as follows: 

p(M(t) =n,M(t + l) = k\Qi(t) = 0, L(t + 1) = L succ ) 

p {Qm(t) > 0,Vm e I,Q r (t) = 0,Vr e I M \I,M{t+l) = k | Q x {t) = 0,L(t+ 1) = L succ ) 

ICIm ■ \I\=n 

£ p (Qm{t) > Q,Vm e I,Q r {t) = Q 7 Vr e I M \I \ Qx(t) = 0,L(t + i) = L succ ) 

ICIm ■ \I\=n 

■ P(M{t + l) = k\ Q m (t) > 0, Vm G /, Q r {t) = 0, Vr e 7 M \ /, Q x {t) = 0, L(t + 1) = L succ ) . 

(C12) 

The second conditional probability inside the summation in the second step in Equation dC. 12b can be expanded 
as follows: (Remember that I C Im and |/| = n.) 

P(M(t + l) - k | Q m (t) > 0,Vm g J,Q r (*) =0,Vr G I M \I,Qi{t) = 0, L(i + 1) = L succ ) 

= P ( Mi - t + !) = fc ' D dt + 1) = 1 | Qm(t) > 0, Vm G I, Q P (t) = 0, Vr G I M \ 7, 
lei 

Q 1 (t) = 0,L(t + l) = L aucc ) 

= j2 P { D l( t + 1 ) = 1 l = ». > > + 1) = ^S"Cc) 

ie/ 

•P(M(t+l) = fc| Qm(i) >0,VmG J ! Q r (t) = 0,VrG/ M \/, 

Qi(t) = 0,i(t + 1) = L MOC ,A(* + 1) = 1) 

= •^(^( t + 1 ) = fc |Qm( t )>o ) v™e/,g r (t) = o,VrG/ M \/, 



Qi(t) = 0, L(t + 1) = L succ , A(i + 1) = 1) 
^ ( I j . p(M(t + i) = k,Qi{t) = 1 | Q m (t) > 0,Vm G J,Q r (t) = 0,Vr G 7 M \/, 



Qi(t) = 0, L(t + 1) = L succ , Dj(t + 1) = 1) 
^ (-) -P(A<(t + l) = *,Qi(t) > 1 | Q m {t) >0,Vm G /,Q r (i) = 0,Vr G J M \/, 



Qi(t) = 0,£(t + l) = £ s „ cc ,A(i + l) = 1).(C13) 



In the first step we condition on the queue from which the departure occurs at t + 1. Note that a departure can 
occur from a queue at t + 1, only if the queue is non-empty at t. Thus, the departure in the success channel slot can 
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occur from one of the queues in /. In the second step we factor out the probability that the departure occurs from the 
(specific) non-empty non-tagged queue at t + 1 given that there are N(t) = n non-empty nodes in the system at t. 
In the third step we apply Equation In the fourth step we condition on whether the non-empty non-tagged queue 
from which the departure occurs at t + 1 contains exactly 1 packet at t. 

The conditional probability inside the first summation in the last step in Equation (IC. 13b can be expanded as 
follows: 

P(M(t+l) = k,Qi(t) = 1 | Q m (t) > 0,Vme/,Q r (t) = 0,Vr 6 I M \I, 

Qi(t) = 0,L(t+ 1) = L SUCCl Di{t + 1) = 1) 
= P(Qi(t) = 1 | Q m (t) > 0,Vm € I,Q r (t) = 0,Vr e I M \I, 

Qx{t) = 0,L(t+ 1) = L succ ,Di{t + 1) = 1) 
• P(M(t + l) = k\ Q m (t) > 0, Vto G / \ {/}, Q,(t) = 1, Q r (t) = 0, Vr e J M \ J, 

Qi(t) = 0,L{t + 1) = L SUCCl Di{t + 1) = 1). (C.14) 

The second conditional probability on the right hand side in Equation (IC. 14b can be simplified as follows: (Re- 
member from Equation (IC. 13b that I, I G I, denotes the index of the non-tagged queue from which the departure 
occurs at t + 1.) 

P(M{t + l)=k\ Q m (t) > 0, Vto € / \ {I}, Qi{t) = 1, Q r (t) = 0, Vr G 7 M \ I, 

Qx(t) = 0, L(t + 1) = L SMCC , + 1) = 1) 
= P(M(t + 1) = * | Q m {t) > 0, Vm G J \ Q,(t) - 1, Q r (t) = 0, Vr G / M \ /, 

Qi(i) = 0, L(f + 1) = L succ , D,(t + 1) = 1, D«(* + 1) = 0, Vu G 7 M \ {?}) 

= P (\ ^2 1 {(Q m (t)+A m (t+l))>0} + 1 {A i (t+l)>0} + ^2 1 {^(t+l)>0} J = k I Qm(t) > 0, 
\me/\{i} rei M \i J 

VtoG J\{/},L(t + l) = L succ ) 

= P( 1{A,(*+1)>0} + X! 1 {^,-(t+l)>0} = - " + 1 | i(< + 1) = Lsucc) 
\ r£l M \I J 

M-n 



k — n + 1 



(1 - s(0)) fc - n+1 s(0) M - fc - 1 . (C.15) 



In the first step we have used the fact that, given that Di(t + 1) = 1, we have D u (t + 1) = 0, Vm G Im \ {I}, 
since at most one departure can occur in a channel slot. In the second step we applied Equations © and (fTJ, and the 
conditioning queue length information of the first step. Note that, Di(t + 1) = 1 and Qi(t) = 1 cancel each other 
in the indicator variable corresponding to /. Thus, the queue from which the departure occurs at t + 1 can be empty 
at t + 1 if it does not receive any packets in the (t + l)*' 1 channel slot. In the third step we have used the fact that 
Q m (t) > 0, Vm G I\{1}, implies that X) me j\{;} 1 {(Q m (t)+A m (t+i))>o} = IA Wl = n-1. The fourth step follows 
by observing that the sum of the M — n indicator variables in the third step is equal to k — n + 1 if and only if arrivals 
occur to exactly k — n + 1 of the corresponding queues, and, given that L(t + 1) = L succ , the probability that arrivals 
occur (resp. do not occur) to a queue is given by 1 — s(0) (resp. s(0)). 

The conditional probability inside the second summation in the last step in Equation (IC. 13b can be expanded as 
follows: 
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P(M(t+l) = k,Qi(t) > 1 | Q m (t) > 0,Vme/,Q r (i) = 0,Vr 6 I M \I, 

Q 1 (t)=0,L(t+l) = L succ ,Di(t + l) = 1) 
= P(Qi{t)>\ \ Q m (t) > Q,Vm e I,Q r (t) = Q,Vr e I M \I, 

Qi(t) = 0, L(i + 1) - i succ , + 1) = 1) 
• P(.M(*+1) = fc | Q m (t) > 0,Vm €/\{/},Qi(*) > l,Qr(*) = 0,Vr 6 7 M \ J, 

Qi(t) = 0, L(f + 1) = i succ , A(i + 1) = 1). (C.16) 



The second conditional probability on the right hand side in Equation ( IC. 16b can be simplified by similar argu- 
ments as follows: 



P(M(t + l) = k\ Q m {t) > 0, Vm G I \ {I}, Qi(t) > 1, Q r (t) = 0, Vr G I M \ P 

Qi(t) = 0,L(t+l) = L auec ,D l (t + l) = 1) 
= P(M(t + l) = k\ Q m (t) >0,VmeI\ {I}, Qi{t) > 1, Q P (t) = 0, Vr e 7 M \ 7, 

Qi(t) = 0,i(f + l) = i succ ,A(t + l) = l,-D u (t + l) = 0,VuG/ M \{0) 



where, in the second step, we observe that (Qi(t) — 1) > 0, and hence, the indicator variable corresponding to I is 
equal to 1 . 

We emphasize that, VI C Im, V/ € I, the conditional probability 



P(Qz(i) = 1 | Q m (t) > 0,Vm G J,Q P (t) - 0,Vr G I M \I,Qi(t) = 0,L(t+l) = L succ , D,(t + 1) = l), 



which appears in Equation dC. 14ft . and the conditional probability 

P(Qz(i) > 1 | Q m {t) > 0, Vm G J, Q r (t) = 0, Vr G I M \ I, Qi(i) - 0,L(i + 1) = L succ ,Dt(t + 1) = l), 

which appears in Equation ( IC. 16b are not known, precisely because the queue lengths of the non-tagged nodes are 
not tracked by the process {X(t)}. Moreover, the above unknown probabilities cannot be obtained from the known 
quantities, namely, the arrival rate A and the state-dependent attempt probabilities /3 n 's. All the other probabilities, 
namely, p id i e ,n, Pcoii,n, Psucc.n, < n < M, d(j), c(j), s(j), j = 0,1,2,.. ., can be obtained from A and the /3„'s. 
Later in this appendix we shall derive the transition probabilities out of a state in which the tagged node is non-empty, 
and it will turn out that the transition probabilities then involve unknown probabilities of the form 



P(Qi(t) = 1 | Q m {t) > 0, Vm G P Q r (t) = 0, Vr G hi \ P Qi(t) =i>0,L(t+l) = L succ , D l (t + 1) = l), 



P(Qi(t) > 1 | Q m {t) > 0,VmG/,Q r (t) = 0,Vr €l M \J,Qi(t) = i> 0,L(t + 1) = L succ , D t (t + 1) = l). 

Recall that I, I G I, denotes the index of the non-tagged queue from which a departure occurs at t + 1, and I 
denotes the set of non-tagged queues that are non-empty at t. At any t, I is a specific subset of Im, and I takes some 




Vm G 7 \ {I}, Qi(t) > l,L(t + 1) = L succ ) 
P ( ^2 1 {A r (t+i)>o} = k - n . | L(t + 1) = L succ ) 




(C.17) 



and 
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specific value between 2 and M. However, due to exchangeability, once we fix the number of non-tagged queues 
that are non-empty at t, i.e., once we fix |/|, the values of the above unknown conditional probabilities do not depend 
on the specific set / and the specific queue I, I £ I, from which the departure occurs at t + 1. Thus, given that 
A4(t) — \I\, we do not have to retain the empty /non-empty status of the non-tagged queues, except that of the queue 
from which the departure occurs at t + 1. We must, retain the fact that Qi (i) > 0; otherwise, Qi (t) = would imply 
Di(t + 1) = 0, which is a contradiction. 

Applying the above exchangeability argument, we can write, Vi > 0, Vz > 0, VI C Im, V/ 6 /, 

P(Qi(t) = 1 | Q m {t) > 0, Vm e /, Q r (t) = 0, Vr 6 I M \ I, Qi(t) = i, L(t + 1) = L succ , D { (t + 1) = l) 
= P{Qt(t) = 1 | Qx{t) = i,M(t) = \I\,Qt(t) > 0, Di(t + 1) = 1, L(t + 1) = L succ ). (C.18) 

To simplify notations, we define, Vt > 0, Vi > 0, Vn, Q < n < M - 1, VZ, 2 < I < M, 

q(i, n, t) := P(Qi(t) = 1 | Q x {t) = i,M(t) = n, Qi(t) > 0, Di(t + 1) = 1, L(t + 1) = L succ ). (C.19) 

Note that, q(i, n, t) on the left hand side of Equation ( IC. 19b does not involve I because, due to exchangeability, the 
value of the conditional probability on the right hand side of Equation (IC. 19b does not depend on the specific value of 
I. Notice that, Vi > 0, Vi > 0, Vn, < n < M - 1, VZ, 2 < I < M, we have 

P(Qi(t) > 1 | Q 1 {t)=i,M{t) = n,Q l {t) > 0,Di(t + l) = l,L(t + 1) = L succ ) 

= l-P(Qi(t) = 1 | Q 1 (t) = 0,M(t)=n,Qi(t)>0,D l (t + l) = 1, L(t + 1) = L succ ) 

= l-q(i,n,t), (C.20) 

Substituting Equation ( IC.15I ) in Equation dC.14t . and applying the notation given by Equation ( IC. 1 9b . we obtain 

P(M(t + l) = k,Qi{t) = 1 | Q m (t) > 0,Vm€l,Q r (t)=0,WeI M \I, 

Qi(t) - 0, Lit + 1) - L succ , D { {t + 1) = 1) 

= q(0,n,t) ■ ( 7 M - n i )(l-s(0)) k - n + 1 s(0) M - k - 1 . (C.21) 



K k - n + 1 

Substituting Equation ( IC. 17b in Equation ( IC. 16b . and applying Equation ( IC.20b . we obtain 

P(M(t + l) = k,Qi(t) > 1 | Q m (t) > 0,Vme/,Q r (i) =0,Vr eJ M \I, 

Qi(t) = 0, Lit + 1) = L succ , D { (t + 1) = 1) 

= {l-q{Q,n,t)) • (^~^~ i yi-s{0)) k - n a(0) M - k - 1 . (C.22) 

Substituting Equations dC.21b and dC.22b in Equation dC.13b . we obtain 

P(M(t + 1) = k | Q m (t) > 0, Vm e /, Q r (t) = 0, Vr 6 I M \ /, Qx(t) = 0, Lit + 1) = L 



\M-k-l 



lei 

^fi) .(l- ff (0,n,t)) • f M ;_ n ji " 1 )(l- S (0)) fc - S (Q) M - fc - 1 



g (0,n,t) • ( L }l n ",)(! s,:u "-'.sue, 

+ (l-g(0 J n,t)) ■ (^""^lO - ".siOi"' ; " '. (C.23) 
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Substituting Equation ( IC.23b in Equation (1C.121 I. we obtain 

P(M(t)=n,M(t + l) = k | Q 1 (t) = 0,L(t + l) = L succ ) 

p {Qm{t) > 0,Vm e I,Q r (t) = 0,W € I M \I | Qi(*) = 0,L(i + l) = L a 

7C/m : |/|=n 

M-n ' 



q(0,n,t) 



Jt-n + 1 
+ (l-q(p,n,t)) 



P(M(t) 



M-n-1 
k — n 

{t)=Q,L{t + l) = L succ ) 
M-n 
k-n + 1 

M-n-1 



(l-s(0)) fe - n+1 s(0) M - fe - 1 

(1 - s(0)) fe -"s(0) M - fe - 1 



9(0, n,t) 
+ (l- ? (0,n,t)) 



(l-s(0))* -B+1 s(°) A '~*~ 1 

(1 - s(0)) fc -"s(0) M - fc - 1 



k — n 



(C.24) 



where in the second step we take the common factor within square brackets outside the summation, and then observe 
that the summation is equal to the conditional probability P(.M(i) = n | Qi(t) = 0,L(t + 1) = L succ }. 
Substituting Equation ( IC.24I ) in Equation ( IC.1U . we obtain 

P(M(t +l) = k\ Q x (t) = 0,M(t) = n, L(t + 1) = L succ ) 
= g(0,n,t). (^^"JCI-^O))*-^ 1 .^-*- 1 



k — n 



+ (l-q(p,n,t)) 

(1 - s(0))' c - n s(0) M -' £ - 1 
+ g(0, n, t) • (1 - s(0))' ! - n s(0) M -' s - 1 ■ 



k-n „fr>\M-k-l 



(i- S (o)r n S (o) 



M-n- 1 

k — n 



M - n 
k-n+1 



(l-,s(0)) 



M -71-1 

k — n 



(C.25) 



*0') 



where in the first step we observe that PL\4(t) — n I Qi(t) = 0, L(t + 1) = L succ ) gets canceled from both the 
numerator and the denominator on the right hand side in Equation ( IC- lib . 
Substituting Equation ( IC.25I ) in Equation ( IC.lOl l we obtain 

P(Qi(t + 1) - j,M(t + l) = k, Lit + 1) = L succ I Q t (t) = 0,X(t) = n) 

^."""^(l-^o))*-*^-*- 1 

k — n J 

+ Psuccn • s(j) ■ q(0,n,t) ■ (1 - S (0)) fe - n s(0) M - ,: - 1 

Substituting Equations ( IC.8I ). ( IC.9b and ( IC.26b in Equation ( IC 1 we obtain 

P(Q 1 (t + l)=j,M(t + l) = k I Qi(t)=0,M(t) = n) 

L ) (Pi<8e,„ • d(j)(l - d(0))' ! - n d(0) M - fe - 1 + PcoII ,„ • c(i)(l - c(0)) fe -"c(0) M -' - 1 

+ p smcc ,„ • - a(o)) fc " B «(o) M " fc_1 ) 



M-n 
fc - n + 1 



;i--(o)) 



M - n - 1 
k — n 



.(C.26) 



M-n - T 
fc — n 



fe— n n/'nl-M — fc— 1 



+ JW.n • • 9(0, n, *)•(!- s(0)) fc -"s(0) 



M-n 
k - +1 



(1- S (0)) 



M-n-1 
k — n 



(C.27) 



which is the transition probability from a state (0, n), < n < M— 1, to a state (j, fc), j = 0, 1, 2, . . ., < fe < M— 1. 
Next, we derive the transition probabilities out of a state in which the tagged node is non-empty. 
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C.2. When the Tagged Queue is Non-Empty in the Initial State 

The transition probability P(Q±(t + 1) = j,M(t + 1) = k \ Qi(t) = i > 0, M(t) = n) corresponds to the 
case when the tagged node is non-empty at the channel slot boundary t and can be expanded according to whether the 
(t + l) th channel slot is an idle, a collision and a success channel slot as follows: 

P(Qi(t+l) =j,M(t+l) = k | Qx(t) = i > 0,M(t)=n) 

= P{Qi(t + 1) = j,M(t + 1) = k, Lit + 1) = L idle | Qt(t) = i > Q,M(t) = n) 
+ P(Q 1 {t + l) =j, M(t + 1) = k,L(t + l) = L coU | Q x {t)=i> 0,M(t) = n) 
+ P(Q 1 (t + l)=j,M(t + l) = k,L{t + l) = L succ | Q x (t) =i>0,M(t) = n). (C.28) 

C.2.1. The Case of Idle Channel Slot 

The first term on the right hand side of Equation (IC.28b can be simplified as follows: 

P(Qi(t + 1) = j,M(t + 1) = k, L(t + 1) = L ldle \Q 1 (t)=i> 0,M{t) = n) 
= P(L(t + 1) = L ldle \Q 1 (t)=i> 0,M(t) = n) 

■ P(A x {t + l)=j- i,M(t + l) = k | Qt{t) = i > 0,M(t) = n, L(t + 1) = L idle ) 
= P(L(t + 1) = Lidie | N(t) = n + 1) • P(M(t + l) = k\ Qi(t)=i> 0,M(t) = n, L(t + 1) = L idle ) 
■P{A x {t + l)=j-i \L(t + 1) = L ldle ) 

d(j - i) ■ ( M 7 n ~\l- d(0)) k -d(0) M - k -\ (C.29) 



n+l ' 



\ k — n 



where we have used the following facts: (a) (L(t + 1) = L id i e ) (D m (t + 1) = 0, 1 < m < M), (b) Q\{t) = i, 
i > 0, and Di(t + 1) = imply that Q x {t + 1) = i + Ax(t + 1), and hence, Qi(t + 1) = if and only if 
+ 1) = j — i, (c) given that L(£ + 1) = Li d i e , the probability that + 1) = j — i, does not depend on 
anything else and is equal to d(j — i), and (d) D m (t + 1) = 0, 2 < m < M, implies that the set of non-empty 
non-tagged queues can only increase, and, given that Ai(t) — n, we have A4(t + 1) = k, if and only if arrivals occur 
to exactly k — n of the M — n — 1 empty non-tagged queues in the idle channel slot, which happens with probability 

(Vn^-dCO))*""^)"-*- 1 . 

C.2.2. The Case of Collision Channel Slot 

The second term on the right hand side of Equation ( IC.28b can be simplified by similar arguments, and we obtain 



Pcoii,n+i-c(j-i)- ( M n l \l-c{Q)) k - n c{Q) M - k -\ (C.30) 



P(Qi(t + 1) = j,M(t + 1) = k, L(t + 1) = i coH | Qi(t) = i > 0,M(t) = n) 

\ — n - 
k — n 

C.2.3. The Case of Success Channel Slot 

The third term on the right hand side of Equation dC.28b can be expanded as follows: 

P(Qi(t+l)=j,M(t + l) = k,L(t + l)=L succ | Q x {t)=i> 0,M(t)=n) 
= P(L(t + 1) = L succ \Qx(t)=i> 0,M(t) = n) 

■ P(Qi{t + 1) = j, M(t + 1) = k | Qi(t) = i > 0,M(t) = n, L(t + 1) = L succ ) 
= fW,„+i • P{Qx(t + 1) = j, M(t + 1) = k | Qi(i) = t > 0, = n, L(i + 1) = i succ ) 
= p S ncc,„+i • P(<9i(t + 1) - j, M(i + 1) = fc, + 1) = 1 | Qi(i) = i > 0,M{t) = n, L(i + 1) = i succ ) 
+ p S n CC ,„+i • P{Qx{t + 1) = j,M(t + 1) = fc,Z»i(t + l) = j Qi(t) = i> 0,M(t) = n, L(t + 1) = i succ ) , 

(C.31) 
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where in the last step we condition on whether the departure occurs from the tagged queue or from a non-tagged 
queue. 

When the departure occurs from the tagged queue, we have 

P(Qi(t + l)=j,M(t+l) = k,Di(t + l) = 1 | Qx(t) = i > 0,M(t) = n,L(t + l) = L succ ) 
= P(D 1 (t + l) = l | Q 1 (t) = i>0,M(t) = n,L(t + l) = L succ ) 

■ P(A 1 (t + l)=j-i + l,M(t + l) = k | Qi{i) = i > 0,M(t) = n,L(t + l) = L succ , D t {t + 1) = l) 
= P(D 1 (t + l) = l | Q 1 (t) = i>0,N(t) = n + l,L(t + l) = L succ ) 

■ P(M(t+l) = k | Qx(t) = i> 0,M(t) = n,L(t + l) = L aucc , D x (t + 1) = l) 
■P(A 1 {t + l)=j-i + l | L(t + 1) = L ldle ) 

' s(j - i + 1) • ( M ~ n ' 1 )(l- S (0)) k - n s(0) M - k -\ (C.32) 



n+1 \ k — n 

where we have used the following facts: (a) given that the tagged queue is non-empty and n non- tagged nodes are 
non-empty at the channel slot boundary t, the departure occurs from the tagged queue with probability ^tj, (b) given 
that Dx(t + 1) = 1 and Qi{t) = i > 0, we have Qi{t + 1) = j if and only if At(t + 1) = j - i + 1, (c) given 
that L(t + 1) = L succ , the probability that Ai(t + 1) = j — i + 1 does not depend on anything else, and is equal to 
s(j + and (d) (Di(t + 1) = 1) => (D m (t + 1) = 0, 2 < m < M), which implies that the set of non-empty 
non-tagged nodes can only increase, and, given that Mt) = n, we have A4(t + 1) = k, if and only if arrivals occur to 
exactly k — n of the M — n — 1 empty non- tagged nodes in the success channel slot, which happens with probability 

C^- 1 )^-m) k - n m M - k - 1 - 

When the departure occurs from a non-tagged queue, we have 

P(Qi(t + l) =j,M(t + l) = k,Di(t + l) = Qi(t)=i> 0,M(t) = n,L(t + l) = L succ ) 

= P{D 1 {t + 1) = | Qx(t) = i > 0,M(t) = n, L(t + 1) = L succ ) 

■ P(A 1 (t + l) = j -i,M(t + l) = k | Q 1 (t)=i> 0,M(t) = n,L(t + l) = L succ ,D 1 (t + l) = 0) 
= (1 - P{D 1 {t + 1) = 1 | Qi(t) = i > 0, N(t) = n + l,L{t + l)= L succ )) 

■ P{M{t + l) = k Qx{t) = i> 0,M(t) = n, L(t + 1) = L succ , D x {t + 1) = 0) 

• P(A X (t + 1) = j - i | L(t + 1) = L ldie ) 

= ■ s(j - i) ■ P(M(t + 1) = k\ Qi(t) = i > 0,JM(f) = n,L(t + 1) = L succ ,D x (t + 1) = 0).(C33) 

?) + 1 v 1 ' 

The conditional probability in the last step in Equation ( IC.33b can be simplified by similar arguments as before, 
and we obtain 

P(M(t + l) = k | Qx{t) =i> 0,M(t) = n,L(t + l) = L aucc ,D 1 (t + l) = 0) 



M-n-l 
k — n 



k-n „(r\\M-k-l 



(1-.(0))*-».(0) 



+ ,(,,,<) ■ (1 - »(0))-»,0)"— ■ (( k M - »(0,) - (";;; ')) , (C.34, 

where notice the similarity with the last step in Equation ( 1C.25I I. However, now we have the unknown probability 
q(i, n, t) instead of q(0, n, t) since the tagged node contains i > packets at t. 
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Substituting Equation ( IC.34b in Equation (1C.331 I. we obtain 

P(Qi(t + 1) = j,M(t + 1) = fc, D x (t + 1) = Qi(t) =i>0 } M(t) = n, L(t + 1) = L succ ) 



n + 1 \ fc — n 



s{3-i)-q{i^t)-{\-si^f- n s{^-^-\\ u _ , , (1-^(0)) 



n + 1 * — n + 1/ \ fc — n 

(C.35) 

Substituting Equations ( IC.32I ) and ( IC.35I ) in Equation ( IC.3U and rearranging, we obtain 

P(Qi (t + 1) = j,M(t + 1) = fc, L(t + 1) = L succ \Qx(t)=i > 0,M(t) = n) 
= ?W,» + i ■ s(j - i) ■ ( M ~_" ~ ^ (1 - S (0)) fe - S (0) M - fe - 1 

+ iWc,»+i (-^—) ■ ( M ~ U ~ X V - aCO))*-"*^)^*- 1 «j - i + 1) - s(j - *)) 



n + 1 / V k — n 
n 

n+1 



fWc,„+i ( — -—r ) ■ s(j - i) ■ q(i, n, t) ■ (1 - s{Q)) k - n s(Q) M ^ 



M-n V , ss /M-n-l\\ 
fc-n+lJ^-^-i fc-n JJ- (C36) 

Substituting Equations ( IC.29I ). dC.30b and dC.36b in Equation ( IC.28b , we obtain 

P(Qi(t + l)=j,M(t + l) = k | Qi(t)=t>0,Al(t) = n) 

L ) (pwie.^+1 • d(j - i)(l - d(0)) fc -"d(0) M - fe - 1 

+PcoH,„+l ■ c(j - - c(0)) fc -"c(0) M ^ fc - 1 

+JW,n+i • «(j - 0(1 - s(0)) fc -"s(0) M - fe - 1 ) 

+ fWc,»+l ( -^-] ■ ( M 7 " ~ \l ~ S (0)) fc -" S (0) M ~ fe - 1 (S(j -Z+l)- S(j - 0) 



M-n - T 
fc — n 



n + 1 J \ k — n 
n 

n+1 



+ fWc,n+i ( — 7— r ) • *(?' - ' ?ft *)•(!- s(0)) fc -"s(0) M - fc - 1 



which is the transition probability from a state (i, n) to a state (j, k), i = 1, 2, . . ., j = 0, 1, 2, . . ., < n, k < M — 1. 
It is important to emphasize again that the transition probabilities of the process {A"(0} involve the unknown 

probabilities q(i, n, i)' s which depend on t. Certain transition probabilities can be zero, since ( ) = when b < 



or a < b, and the probabilities d(j), c(j) and s(j) are zero for j < 0. 



D. Transition Probabilities of the Process {X(t)} 



The transition probabilities of the process {X(t)} can be derived from the transition probabilities of the process 
{X (t)} as follows. First we apply Lemma |4~T1 and redefine q(i,n,t), Vt > 0, Mi > 0, Vn, < n < M — 1, Vi, 
2 < I < M, as given by Equation ( TTOb . Then, we apply Approximation l4~Tl and Equation ( fT2b to write 



q(i,n,t) 



q(n, t) if i = 0, 

g(n + l,0 if « > 0. 
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Finally, we replace the q(n, t)'s by the constant time-independent probabilities q(n)'s and obtain the transition prob- 
abilities of the process {X(t)} in terms of the q(n)'s as follows. 

P(Q 1 (t + l)=j,M(t + l) = k | Q x (t) =0,M(t)=n) 

= ( M ~™ ~ ^ (Pidie,n ■ d(j)(l - d(0)) fc -"d(0) M - fc - 1 + Vcoll>n ■ c(j)(l - c(0)) fe -"c(0) M - fc - 1 

+ Vsucn ■ S(j)(l - s(0)) fc -" S (0) M - fc - 1 ) 

+ Psu CC , n ■ s{j) ■ ?(n) • (1 - «(0))*-" fl (0)^-*- 1 • (( fc M ~ - S (0)) - ( M ~" ~ *)) , (D.l) 

P(Qi(* + l) = j,.M(t + l) = fc I Qi(i) = • > 0,M(t) = n) 

= ( M fc ~_ n ~ fr^-^ 1 ■ - W 1 - rf(o))^ n d(o) M -' £ - 1 

+p GO H,„+i • c(j - i)(l - c(0)) fc -"c(0) M - fe - 1 

+2Wc,„+i • s(j - i)(l - S (0)) fe -" S (0) M " fe - 1 ) 

+ Psuccn+l (-1—) ■ H ~ 1 )(l- S (0)) fe -" S (0) M - fe - 1 (S(j - i + 1) - S (j - 0) 



rt + 1/ V k — n 



M-k-1 



M — n \ , . {M — n — l\\ 

Due to the constant time-independent probabilities q(n)'s, the process {X(t)} is a time-homogeneous DTMC 
embedded at the channel slot boundaries. However, the constant time-independent probabilities q(n)'s are not known, 
and they must be treated as unknown parameters of the process {X(t)}. 

E. Derivation of Equation (f30t 

Recall that the DTMC {Q^ (t),t > 0) i s stationary and ergodic. We apply the discrete time version of Birkoff's 
strong ergodic theorem (see page 800 of 12111 ) to obtain p( d > (j), < j < K — 2. (Note that, the time argument t in 
Theorem lE. H is discrete.) 



Theorem E.l (Birkoff's Strong Ergodic Theorem). Let {X(t),t > 0}, be a stationary and ergodic process, and let 
/(•) be a function that maps realizations of the process (i.e., {X(t, u>),t > 0}) to the set of real numbers K, such that 
E(\f({X(t),t > 0})|) < oo. Then, with probability 1, 

1 r_1 

lim - V f({X(t + 1% t>0}) = E(f({X(t), t > 0})) . ■ 

r— voo T — * 

t'=0 

Noting that departures can occur only at the end of channel slots, we write 

T-1 t-1 

1 {D< K '(to+l+t')=l} 1 {oi K) (to+l+t')=j} 1 {z5< fc) (t +i+t')=i,Qr ) (to+i+t')=i} 
p^(j) = lim 4 -° ' t,=0 



E 1 {^( to +l+t')=l} ™ t ^ 1 {^ ) (*o 

*'— n J + '— n ^ 



t-1 

+ l+t') = l} 



f'=0 

where, applying ergodicity of the process {Q^(0,i > 0}, we equate the probability p^ d \j) with the long-term 
fraction of departures from the tagged queue that leave j packets behind (starting from an arbitrary time to + 1). Then, 
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we can apply Theorem lE.il observing that expectation of the indicator function of an event is equal to the probability 
of the event, and hence, is always bounded between and 1 . Thus, we have 



E ( 1 {D[ K \t ()+ i)=i,Q[-\t 0+ i)= 3 } ) _ p( D [ K \t Q + 1) = l,Q{ K) (t + 1) = j) 



E( 1 {^(t +l)=l} 



P(D[ K \t + l) = l) 



K M-l 

E E P (Q{ K) (to) = i,M^(t ) = n, D[ K \t + 1) = 1, Q{ K \t + 1) = j) 

i=l n=0 

K M-l 

E E HQi K) Qo) = ^M^ K \t,) = n,D[ K \t + !) = !) 

i=\ 71 = 



(E.l) 



where, in the last step, we expand by summing over all possible states at time to- Notice that, in the last step, we do 
not consider the possibility i = since (^D[ K \t a + 1) = lj \Q{ K \t ) > o). 
For < j < K — 2, Equation (IE. U can be simplified as follows: 



k M-i j P{Q{ K) (t ) =i,MW(to) =n) \ 
EE x p i D i K) ^ + !) = 1 I Q^fo) = *,-M W (to) - n) 
<= 1 »=° V x P(Q^Vo + l)=j | Q[ K) (to) =i s A^W(*o) = n s I>i K) (to + 1) = l) / 



K M-l 



P(Q[ K \t )=i,MW(t ) = n) 



i=l n=0 

K 



\ x P(D[ K \t + 1) = 1 | Q[ K \t ) = i,M^ K \t ) = n) 



m-i ( P{Q{ K) (t ) =i,M {K) {to) =n) \ 
EE X P(D[ K \to + 1) = 1 I Q{ K \to) = hMWfa) = n) 

1=1 V x P(A[ K) (t a + l)=j-t + l\Q[ K \t )= i ,M( K Hto)=n,D[ K \t + l) = l) J 



k ' < P(Q[ K \t ) = i,M^(t ) = n) 

x P(D[ K \t + 1) = 1 | Q[ K) (t ) = i,M {K) {t a ) = n) 



i=l n =Q 



EE i 

K M-l 

E E P(Qi K) (t°) = ^ W (to) = n) ( P -^f ) s(j 

i=l n=0 J- 

K M-l , 

E E p (Q ( i K ) (*o) = ^^ w (t ) - ^ (>'-■ 

i+iM-i / 

EE- W (^)(^f ) S (i-i+D 



n + 1 



= 1 71=0 



if M-l , . 



(E.2) 



i = l 71 = 



In the second step, we used the fact that, given Q\ K ^ (to) = i and (to + 1) = 1, we have Q± (to + 1) = j 
if and only if A[ K \t + 1) = j — i + 1. In the third step, we observe that Qi (to) = i > and A4 (if) (t ) = n 
implies that N^ K ^(to) = n + 1. Thus, the probability that the (to + l) th channel slot is a success channel slot 
is given by p SU cc,n+i, and given that a success occurs, the departure occurs from the tagged node with probability 
l/(n + 1). In the fourth step, we let to — > oo on both sides. Since, the left hand side of Equation ( IE.2I ) does not 
depend on to, it remains unchanged. However, for the stationary and ergodic DTMC {Q^ K ' (t),t > 0}, we have, 
lim t0 _ i . oo P(Q[ K \to) — i, A4( K ^(to) = nj = ir( K ^(i,n) where ir^ K \i,n) denotes the stationary probability that 
the tagged queue contains i packets and n non-tagged queues are non-empty. Notice that the first sum index of the 
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numerator in the last step on the right hand side of Equation ( IE. 2b runs only up to j + 1. For i > j + 1, the argument 
j — i + 1 of s(j — i + 1) becomes negative so that s(j — i + 1) = 0. 
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